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I.  INTRODUCTION 


Directional  hydrophone  arrays  are  often  used  to  determine  the  physical 
attributes  of  a  distant  source  or  reflector  of  acoustic  waves.  Although  the 
range,  velocity,  and  bearing  of  the  target  are  often  of  greatest  Interest, 

In  many  cases  It  Is  Important  to  learn  as  much  as  possible  about  Its  geometric 
structure,  and  to  determine  Its  precise  bearing  angle  and  depression- 
elevation  angle  with  respect  to  background  objects.  A  single  beamformer 
may  be  used  for  this  purpose,  but  better  results  can  often  be  obtained  by 
using  two  beamformers,  for  example,  one  having  a  sharp  peak  and  the  other 
a  null  at  broadside. 

Complementary  beams  have  long  been  known  to  be  useful  for  precision 
acoustic  or  radio  direction  finding.1  Although  more  sophisticated  tech- 
nlques  have  been  developed  In  recent  years,  *  2-beam  techniques  still  have 
value.  They  provide  a  resource  that  has  not  yet  been  exhausted,  particularly 
with  regards  to  broadband  signals,  and  they  are  much  more  amenable  to 
performance  prediction  than  other  techniques. 

In  this  paper,  some  familiar  2-beam  techniques  are  extended  to  the 
broadband  case.  In  particular,  the  amplitude  comparison  concept  Is 
extended  to  a  scale  comparison  method  that  Is  applicable  to  signals  of 
arbitrary  spectra  and  angles  of  Incidence.  Moreover,  several  new  results 
are  given  that  are  of  Interest  even  In  the  single  frequency  case.  In 
contrast  to  most  recent  studies,  two  outputs  are  produced  that  are  related 
by  a  simple  proportionality  rule  rather  than  by  a  time  delay.  This 
approach  makes  It  possible  to  use  simpler  statistical  models  for  the 
corrupting  effects  of  multipath  distortion  or  complexities  of  target 
structure.  An  appendix  Is  Included  on  phase  front  distortion  as  It  affects 
estimates  of  wave  source  direction. 
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The  report  Is  written  In  tutorial  style,  but  It  contains  several 
new  results  and  methods.  We  begin  with  a  review  of  some  basic  properties 
of  arrays  as  seen  from  a  signal  processing  point  of  view. 
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II.  SIGNAL  PROCESSING  CHARACTERISTICS 
OF  A  LINE  ARRAY  BEAMFORMER 


Suppose  a  plane  wave,  not  necessarily  sinusoidal,  arrives  at  angle 
0  as  shown  in  Fig.  1.  Note  that  e  is  the  angle  between  the  wavefront 
and  the  positive  x-axis  as  well  as  the  angle  by  which  the  wave  source  is 
displaced  from  broadside  aspect.  If  the  instantaneous  acoustic  pressure 
at  the  origin  is  p(t),  then  at  an  arbitrary  point  along  the  x-axis.  Its 
measured  value  must  be  p[t-(x/c)sine],  where  c  Is  the  speed  of  propaga¬ 
tion.  Suppose  a  receiving  aperture  with  aperture  shading  function  w(x) 
Is  ured  to  sense  the  wave  field  along  the  x-axis  so  as  to  produce  an 
output 


q(t)  ■/ w(x)  p(t 


(1) 


(Integrals  with  unlabeled  limits  will  be  assumed  to  extend  over  the  range 
-oo  to  +«.)  The  aperture  shading  function  w(x)  may  be  effected  by  the 
physical  characteristics  of  the  receiving  transducers,  by  electronic  cir¬ 
cuitry,  by  digital  computation,  or  by  a  combination  of  those  methods. 

In  any  case,  the  entitles  that  perform  the  spatial  Integration  of  Eq.  (8) 
will  be  collectively  referred  to  as  a  beamformer.  Used  here,  this  term 
does  not  Include  beamsteering,  a  processing  task  that  will  be  neglected 
in  our  treatment. 

For  an  aperture  of  finite  length,  the  function  w(x)  is  zero  outside 
some  finite  region,  and  the  integration  limits  can  be  collapsed  to  that 
domain  if  desired.  The  shading  function  w(x)  is  permitted  to  take  nega¬ 
tive  values,  but  the  integral  of  its  magnitude  must  be  finite.  In  the 
case  of  a  receiving  array  composed  of  a  finite  number  of  point 
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transducers,  w(x)  Is  just  a  series  of  Impulses.  If  we  substitute  kt  for 
t  In  Eq.  (1),  with  k  defined  as  cosine,  the  result  Is 

q(Kt)  =  /w(x)  p[k(t-x)]  dx  =  w(t)  *  p(kt)  ,  (2) 

where  *  denotes  convolution.  But  convolutions  obey  the  following  scaling 
property,  as  can  easily  be  verified:  if  f (t)=g(t)*h(t) ,  then 

f (at)  =  j a {  g(at)  *  h(at)  .  (3) 

Application  of  this  property  to  Eq.  (2)  with  1/k  as  the  scalar  a  and 
q(kt)  as  f(t)  gives  the  result 

q(t)  =  w0(t)  *  p(t)  ,  (4) 

where  the  temporal  aperture  function  for  direction  e  is  defined  as 

we(t)  =  l<|  w(l<)  =  Irniel  w(slW)t  '  ^ 

This  new  function  may  be  visualized  as  an  area-preserving  projection  of 
w(x)  onto  a  t-axis  that  is  directed  along  the  path  of  wave  propagation 
(see  Fig.  2).  The  Fourier  transforms  of  w(t)  and  w0(t)  must  then  obey 
the  relationship  W. (f)=W[(f/c)s1ne].  Furthermore,  W.(f)  has  the  constant 
value  Mq  when  es0,  where  MQ  is  the  area  under  the  curve  described  by  w(x), 
i.e.,  it  is  the  "zeroth"  moment  integral  of  the  aperture  shading  function. 
Equation  (4)  gives  the  time  domain  representation  of  the  signal  pro¬ 
cessing  function  of  the  beamformer.  Taking  the  Fourier  transform  gives 

Q(f)  -  WQ(f)  P(f)  =  w(£  sine)  P(f)  ,  (6) 
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FIGURE  2 

TEMPORAL  APERTURE  FUNCTION 


which  Is  the  frequency  domain  representation  of  the  beamformer's  signal 
processing  behavior.  On  the  other  hand,  If  Eq.  (1)  Is  attacked  by 
expanding  p[t-(x/c)s1ne]  In  a  power  series,  i.e., 


p(t  -  |  sine)  =  p(t)  -  sine  £  p(t)  +  ...  ,  (7) 


then  the  result  is 

q(t)  =  M0p(t)  -  sine(M^/c)  p(t)  +  ...  ,  (8) 

where  the  nth  term  in  the  series  is 

M 

-  -I}-  (sine)n  c“n  p(n)(t)  ,  (9) 

n  • 

with  Mn  being  the  nth  moment  integral  of  the  aperture. 


M 


n 


dx 


(10) 


Note  that  the  definition  of  MQ  is  consistent  with  that  given  previously, 
and  W( 0 )=Mq .  Equation  (8)  gives  the  derivative  series  representation 
for  the  signal  processing  behavior  of  the  beamformer.  An  important  con¬ 
sequence  of  this  representation  is  that  as  9  approaches  zero  (broadside 
aspect)  the  output  becomes  like  M0p(t),  which  means  that  if  MQ/0  then  the 
dominant  term  has  no  first  order  dependence  upon  6.  On  the  other  hand, 

If  Mq=0  and  M^Othen  the  dominant  term  in  the  output  is 
(M-|/c)p(t)  sine  ~  (-M^/c)p(t)e .  The  resulting  first  order  dependence  up¬ 
on  0  of  the  dominant  term  makes  this  type  of  beamformer  useful  for  esti¬ 
mating  the  direction  of  plane  waves  arriving  at  near-broadside  aspect, 
if  the  derivative-taking  action  of  the  array  is  taken  into  account. 
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To  emphasize  the  directional  characteristic 
particular  frequency  f  ,  It  is  useful  to  rewrite 


of  the  output  at  a 
Eq.  (6)  in  the  form 


Q(f0)  -  K(e.f0)  P(f0) 


(11) 


where  K(e »fQ)  is  the  directivity  function  for  frequency  f  , 

K(e,f0)  £  W[(f0/c)sine]  .  (12) 


This  complex  function  relates  the  amplitude  and  phase  of  the  output  com¬ 
ponent  at  frequency  f  to  that  of  the  incident  pressure  wave  p(t).  In 
particular,  if  p(t)  is  a  pure  sinusoid,  i.e., 


p(t) 


j2trf  t 
Pe  0 


» 


where  P  is  the  complex  phase  amplitude  factor,  then 


P(f)  =  l/2[Pfi(f-f0)  +  P*6(f+f0)] 


(13) 


(14) 


which,  after  substituting  into  Eq.  (6)  and  taking  inverse  transforms, 
leads  to 


q(t)  =  Re 


Qe 


j  2irf  Qt 


(15) 


where 

Q  4  K(e,f0)P. 


(16) 


Equation  (16),  used  here  to  define  the  beamformer  output  phaser  Q, 
might  also  be  regarded  as  a  working  definition  of  the  directivity  func¬ 
tion  K(e,fQ),  given  in  terms  of  the  complex  envelope  of  the  beamformer 
output  for  an  incident  sinusoidal  wave.  The  directivity  function  is  often 
expressed  with  its  broadside  value  normalized,  i.e.,  as  K(6,f0)/M0,  in 
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the  literature.  But  this  cannot  be  done  If  Mo°0,  l.e..  If  the  area 
under  w(x)  Is  zero,  so  that  the  directivity  pattern  has  a  null  at  broad¬ 
side,  and  therefore  no  normalization  will  be  applied  here. 

It  Is  worth  pointing  out  that  Eq.  (12)  underpins  the  following  well- 
known  Interpretation  of  spectral  transforms:  any  graphical  compendium 
of  Fourier  transforms  may  be  regarded  as  a  compendium  of  directivity 
functions  If  one  merely  marks  the  frequency  axis  according  to  the  non¬ 
linear  angular  scale  e=arcsin(cf/f0) ,  where  f  Is  the  frequency  variable 
In  the  graphical  representation  of  the  transform.  One  needs  only  that 
part  of  the  plot  for  which  |cf/f  [<1,  i.e.,  |f|<f0/c. 
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III.  SYMMETRY  PROPERTIES 

Suppose  an  aperture  shading  function  w(x)  has  even  symmetry, 
w(x)sw{-x).  Then  Its  odd  moment  Integrals  M-j.M^,...,  vanish,  which 
Implies  via  the  derivative  series  representation  of  Eq.  (8)  that  the 
beamformer  output  must  be  composed  of  p(t)  and  Its  even  derivatives, 
with  only  even  powers  of  9.  The  Fourier  transform  of  w(x)  Is  iecessarlly 
real-even,  which  means  the  directivity  function  K(e,fQ)  is  real  valued 
and  possesses  even  symmetry  with  respect  to  both  of  the  variables  e  and 
f  .  The  beamformer  output  carries  no  information  as  to  the  sign  of  e, 
and  all  of  its  sinusoidal  components  are  in  phase  (allowing  for  reversed 
sign)  with  those  of  the  acoustic  pressure  at  the  origin.  At  broadside 
aspect  the  output  Is  relatively  insensitive  to  small  variations  In  the 
angle  of  the  incidence,  exhibiting  only  a  second  order  dependence  on  e. 

By  similar  reasoning,  a  shading  function  with  odd  symmetry, 
w(x)=-w(-x) ,  has  M0=M2=,. . .,=0,  so  its  output  q(t)  is  composed  of  only 
odd  derivatives  of  p(t)  with  odd  powers  of  0,  and  the  Fourier  transform 
of  w(x)  Is  imaginary-odd.  The  directivity  function  K(e,f0)  displays  odd 
symmetry  with  respect  to  0  and  f  .  The  beamformer  output  reverses  In  sign 
whenever  the  direction  angle  0  crosses  broadside  aspect,  and  at  broadside 
displays  a  first  order  sensitivity  to  9  that  is  determined  by  the  size  of 
the  first  moment  Integral  M^ .  Every  sinusoidal  component  of  the  beamformer 
output  q(t)  Is  90°  out  of  phase  with  its  antecedent  component  In  the 
incident  wave  pressure  p(t),  with  the  direction  of  that  phase  shift 
depending  upon  the  sign  of  0. 

Since  the  odd  or  even  symmetry  of  the  aperture  shading  function 
implies  a  similar  symmetry  in  the  directivity  function,  the  terms 
"odd  beamformer"  or  "even  beamformer"  are  entirely  appropriate. 
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IV.  TEMPORAL  INTEGRATION  OF  BEAMFORMER  OUTPUT 


Equations  (4)  and  (5)  can  be  combined  to  express  the  beamformer 
output  more  directly  for  a  general  aperture  shading  w(x), 

q(t)  =  | k| “1  w(t/k)  *  p(t)  ,  (17) 

where  kcc’^sin  as  before.  Now  suppose  q(t)  is  integrated  from  -»  to 
t  by  an  operational  amplifier  circuit.  This  temporal  integration  can  be 
symbolically  represented  as  u(t)*q(t),  where  u(t)  is  the  unit  step 
function.  But 


u(t)  *  q(t)  *  |k|-1  u(t)  *  w(t/k)  *  p(t) 

=  [kj'1  f(t)  *  p(t),  (18) 

where  f(t)  =  u(t)*w(t/k).  Moreover,  the  convolution  scaling  property  of 
Eq.  (3)  implies  that  f(kt)s|k|u(kt)*w(t) .  If  we  define  a  new  aperture 
shading  function 


w.(x)  4  sgn(k)  u(kx)  *  w(x)  ,  (19) 

V 

then  f(t)«kw  (t/k)  which,  when  substituted  into  Eq.  (18),  gives 

V 

u(t)  *  q(t)  =  k  |k|_1  wc( t/k)  *  p(t)  .  (20) 

Except  for  the  multiplicative  factor  k,  the  right-hand  side  of  this 
equation  has  the  same  form  as  that  of  Eq.  (17),  which  means  we  can  rewrite 
it  in  the  equivalent  form  of  Eq.  (1), 
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The  upshot  Is  that  whenever  the  output  of  a  beamformer  having  aperture 

shading  w(x)  Is  passed  through  an  Integrating  circuit,  the  resulting  signal 

Is  Identical  to  what  would  have  been  produced  by  the  companion  beamformer 

whose  aperture  shading  w„(x)  Is  defined  In  Eq.  (19),  except  that  It  Is 

multiplied  by  k“c  sine.  Although  the  definition  of  w.(x)  would  seem  to 

C  -1 

depend  on  the  angle  of  arrival  through  Its  dependence  upon  k°c  sine,  It 
really  depends  only  on  the  sign  of  e,  since  u(kx)=u(x)  If  k>0,  and 
u(kx)-l-u(x)  If  k<0  .  It  follows  that 


t  s1"«) 


dx 


(21) 


u(.t)  *  q(.t)  *  k 


/  w< 


(x)  p 


If  e>0,  but 


wc(x) 


wc(x) 


w(x)  dx 


+00 

w(x)  dx 


(22) 


(23) 


If  e<0. 

If  w(x)  has  odd  symmetry,  then  these  two  formulas  are  completely 
equivalent.  Thus  to  any  odd  beamformer  one  can  associate  a  unique,  even, 
companion  beamformer  whose  aperture  shading  function  Is  just  the  Integral 
of  the  odd-symmetric  function  w(x)  from  -»  to  x.  Odd  and  even  aperture 
shading  functions  having  this  relationship  are  said  to  satisfy  the  "Kerr- 
Murdock  conditions"  In  the  antenna  literature. For  a  given  odd- 
symmetric  shading  function,  the  corresponding  even- symmetric  shading 
function  so  formed  will  be  termed  the  K-M  companion,  denoted  WKM^' 

Equation  (21)  makes  It  possible  to  determine  the  performance  of  an 
odd  beamformer  from  that  of  a  hypothetical  even  beamformer.  The  directional 
characteristics  of  an  even-symmetric  aperture  are  much  more  familiar  to 
array  designers,  and  to  make  the  transition  to  an  odd  beamformer  It  Is 
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only  necessary  to  Integrate  Its  aperture  shading  to  get  the  even- 
symmetric  aperture  shading  of  Its  K-M  companion  even  beamformer,  whose 
directional  characteristics  will  exactly  predict  those  of  the  odd  beam- 
former  except  for  a  temporal  Integration  (which  applies  a  90°  phase 
shift  to  each  sinusoidal  wave  component)  and  a  multiplicative  factor 
k»c"*s1n6  (which  gives  the  beamformer  Its  characteristic  odd- symmetric 
directivity  function). 

A  more  Important  application  of  this  result,  however,  lies  In  the 
design  of  beamformers  for  precision  direction  finding. 


V.  DIRECTION  FINDING  BY  TRANSFORMS 


The  Fourier  transform  relationship  between  the  beamformer  output 
q(t)  and  the  Incident  wave,  represented  by  the  acoustic  pressure  at  the 
origin,  Is  given  by  Eq.  (6),  repeated  here  for  convenience: 


Q(f)  -  P(f)  w(£  sine)  .  (24) 

Since  q(t)  Is  measured,  c  Is  a  physical  constant,  and  w(x)  is  a  designer's 
choice,  the  only  unknowns  Implicit  In  Eq.  (24)  are  the  Incident  waveform 
p(t)  and  Its  angle  e.  For  the  rest  of  this  report  we  shall  be  concerned 
with  determining  the  Incident  wave  angle  e. 

If  p(t)  could  somehow  be  divined,  then  the  three  transforms  Q(f), 
W(f),  and  P(f)  would  need  to  be  computed  only  once,  and  e  could  be  found 
to  achieve  equality  in  Eq.  (24).  But  In  practice  p(t)  Is  never  known 
beforehand.  Even  if  It  is  produced  as  a  reflection  of  a  known  projected 
wave  (l.e.,  by  active  sonar)  or  by  a  distant  transmitter  of  known  output, 
then  at  best  one  knows  p(t)=ApQ(t-a) ,  where  P0(t)  Is  known  but  a  and  A 
are  unknown.  Recalling  that  the  magnitude  of  a  Fourier  transform  Is 
Invariant  with  respect  to  a  time  shift,  we  may  conclude  from  Eq,  (24) 
that 


I  Q(f )  |  =  I A  (  |PQ(f)|  |w(£s1ne)  |  ,  (25) 

which  may  be  amenable  to  solving  for  e.  Since  there  are  two  unknowns, 

| A |  and  the  Incident  wave  angle  e,  It  is  essential  to  treat  at  least  two 
frequencies.  Unfortunately,  the  solution  for  o  will  be  ambiguous  as  to 


sign,  since  |W(f)|  Is  an  even  function  when  w(x)  is  real,  but  the 
Information  gained  about  e  may  be  useful  despite  that  shortcoming. 

The  "phase- Ignorant"  method  described  above  may  have  some 
applications,  but  It  Is  far  better  to  use  two  beamformers  having  differ¬ 
ent  aperture  shading  functions,  w-|(x)  and  w2(x),  that  spatially  Integrate 
the  acoustic  field  simultaneously  along  the  same  measurement  axis.  In 
practice,  this  may  be  achieved  by  using  two  sets  of  shading  resistors  for 
the  same  array  of  transducers,  by  placing  two  arrays  end  to  end  (or  even 
side  by  side  as  an  approximation),  splitting  an  existing  array  at  the 
middle,  processing  the  individual  transducer  outputs  in  a  computer,  or 
by  a  variety  of  other  means.  Thus  two  outputs,  q-|(t)  and  q2(t),  can  be 
produced,  giving  two  transform  equations  of  the  form  of  Eq.  (24),  from 
which  the  common  term  P(f)  can  be  eliminated,  leaving 

Ql  (f)  W2(£  sine)  =  Q2(f)  «,(£  sine)  ,  (26) 

which  may  be  solved  for  e. 

This  computational  transform  method  for  direction  finding  can  be 
summarized  as  follows.  The  two  beamformer  outputs  q-j(t)  and  q2(t)  are 
Fourier  transformed  to  get  Q^(f)  and  Q2(f),  which  may  be  evaluated  at 
some  set  of  discrete  frequencies  (hence,  a  fast  Fourier  transform  is 
acceptable).  In  order  to  test  Eq.  (26)  for  equality  at  different  values 
of  e,  the  Fourier  transforms,  W^(f)  and  W2(f),  of  the  shading  functions 
must  be  computed  with  (f/c)sin9  substituted  for  f.  This  is  most  easily 
done  if  Wj(f)  and  (f )  have  relatively  simple  functional  forms  comprising 

algebraic  and/or  elementary  transcendental  functions.  Then  e  can  be 
adjusted  to  minimize  the  error  of  Eq.  (26)  over  an  appropriate  range  of 
frequencies  by  any  of  various  criteria,  such  as  minimizing  the  magnitude- 
squared  equation  error.  The  use  of  just  a  single  frequency  may  be 
sufficient. 
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Although  this  method  seems  straightforward  enough »  Its  success  will 
hinge  upon  such  factors  as  the  dissimilarity  of  the  two  beamformers,  the 
distribution  of  signal  energy  In  the  frequency  spectrum,  and  the  effects 
of  the  windowing  that  are  Inevitable  when  computing  the  Fourier  trans¬ 
forms  of  the  beamformer  outputs.  Just  how  these  factors  will  affect  the 
results  Is  difficult  to  determine. 

There  are  ways,  however,  to  extract  e  directly  from  a  pair  of 
beamformer  outputs.  These  methods  are  highly  structured  and  have  a  more 
easily  predicted  performance  than  the  computational  transform  method. 
They  rely  upon  the  use  of  specific  choices  for  the  aperture  shading 
functions  w^(x)  and  w2(x),  or  upon  restrictive  assumptions  concerning 
p(t)  (e.g.,  that  It  Is  sinusoidal).  The  balance  of  this  report  will  be 
devoted  to  a  study  of  these  methods.  However,  one  should  keep  In  mind 
that,  after  performance  predictions  have  been  made,  there  Is  always  the 
option  of  analyzing  the  beamformer  output  using  the  computational  trans¬ 
form  method  just  described,  and  perhaps  of  capitalizing  upon  the  beam- 
former  structure  to  improve  the  computational  accuracy  or  to  compensate 
for  known  Imperfections  in  beamforming. 
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VI.  DIRECTION  FINDING  BY  TIME  COMPARISON 

If  the  aperture  weighting  functions  of  a  pair  of  beamformers  are 
related  by  an  axial  shift,  w^xJ^fx-A),  then  the  beamformer  outputs 
exhibit  a  differential  time  shift,  q^(t)aq2(t-As1ne/c),  as  may  be  verified 
by  expressing  q^ (t)  In  the  Integral  form  of  Eq.  (1)  and  changing  the 
variable  of  integration  to  x'®x-A.  Such  beamformers  may  be  termed  time 
comparison  beamformers  since  the  direction  angle  e  of  the  Incident  plane 
wave  reveals  Itself  as  e^arcsinfcTA"1 ) ,  where  T  is  the  relative  time 
shift.  Axially  shifted  aperture  weighting  functions  may  be  realized 
quite  simply,  for  example,  by  splitting  a  transducer  array  at  the  center 
or  by  placing  two  arrays  end  to  end.  Procedures  for  estimating  the 
relative  delay  between  two  signals,  even  in  the  presence  of  noise,  have 
been  thoroughly  treated  In  the  literature  (a  special  Issue  of  the  IEEE 
Transactions  on  Acoustics,  Speech,  and  Signal  Processing  was  devoted  to 
that  subject).7  A  statistical  crosscorrelation  of  the  two  beamformer 
outputs  may  be  all  that  Is  needed,  at  least  In  theory.  The  problem  Is 
not  an  easy  one,  however,  particularly  If  the  incident  wave  angle  e 
varies  so  rapidly  that  the  analysis  must  be  confined  to  a  brief  sample 
of  output  data.  Although  the  relationship  between  the  beamformer  out¬ 
puts  Is  linear,  it  is  not  memoryless.  Moreover,  If  the  incident  wave  Is 
periodic  with  period  T  ,  then  T  can  be  determined  only  to  within  an 
unknown  Integer  multiple  of  T  ,  resulting  in  ambiguous  solutions  for  e, 
unless  A<Tqc. 

Since  this  ground  has  been  explored  extensively  by  many  Investigators, 
we  shall  confine  our  attention  to  pairs  of  beamformers  that  exhibit  a 
linear,  memory less  relationship  between  their  outputs. 


VII.  DIRECTION  FINDING  BY  SCALE  COMPARISON 


Consider  an  odd-even  pair  of  beamformers  whose  aperture  shading 
functions  satisfy  the  Kerr-Murdock  conditions,  where  w^(x)  Is  the  odd- 
symmetric  aperture  shading  function  and 


wkm(x) 


x 

w-j(x)  dx 


(27! 


As  previously  shown,  the  temporally  Integrated  output  of  the  odd  beam- 
former. 


q^t)  £  u(t) 


-00 


(28) 


Is  related  to  the  unintegrated  output  of  the  even  (K-M  companion)  beam- 
former  by 


?1  (t)  =  kq^t)  a  c'1  sine  qKM(t)  .  (29) 

Thus,  the  incident  wave  angle  e  can  be  determined  unambiguously  over  the 
range  -tt/2  to  +tt/2  by  finding  the  factor  k  that  relates  the  beamformer 
outputs  and  then  setting  e=arcsin(kc) .  Interestingly,  if  the  wave  Is 
presumed  to  have  come  from  some  remote  point  at  range  d,  then  the  dis¬ 
placement  of  that  point  from  the  broadside  plane  (i.e.,  the  crossrange) 
Is  just  kcd,  as  follows  from  elementary  trigonometry.  Determining  the 
value  of  k  can  be  as  simple  as  dividing  ( t)  by  q^t),  provided  that: 
(1)  the  slgnal-to-noise  ratio  (S/N)  is  very  large, 


(2)  the  Incident  wave  Is  truly  planar,  and 

(3)  the  beamformers  have  been  Implemented  with  no  compromises  In 
aperture  shading. 

The  use  of  an  X-Y  oscilloscope  display  Is  particularly  appropriate  for 
determining  the  scaling  factor  that  links  the  two  signals. 

Figure  3  gives  several  possible  pairs  of  aperture  shading  functions 
that  satisfy  the  Kerr-Murdock  conditions,  but  in  practice  they  may  be 
difficult  to  realize  precisely.  The  easiest  way  to  create  an  odd-even 
beamformer  pair  with  outputs  q-j (t )  and  q2(t)  is  to  choose  an  even  aper¬ 
ture  shading  function,  split  the  aperture  at  the  center  to  form  separate 
left  and  right  outputs,  and  then  take  q^(t)  as  the  difference  of  those 
outputs  and  q2(t)  as  the  sum.  In  that  case,  w2(x)s±sgn(x)w1 (x) .  Indeed, 
any  pair  of  beamformers  designed  for  the  time  comparison  technique  dis¬ 
cussed  previously  can  be  converted  to  such  an  odd-even  pair  by  this  sum- 
and-dlfference  method.  Unfortunately  the  only  odd-even  pair  of  shading 
functions  that  satisfies  both  the  Kerr-Murdock  condition  and  the  split 
aperture  condition,  i.e.,  for  which  w1  (x^tsgnfxjw^'fx) ,  is  the  exponen¬ 
tial  pair  given  In  Fig.  3,  which  requires  a  line  array  of  infinite  length 
(admittedly,  the  exponential  aperture  might  be  approximated  by  one  of 
finite  length).  Even  if  one  is  willing  to  forego  the  simplicity  of  a 
split  aperture  beamformer,  It  must  be  appreciated  that  the  Kerr-Murdock 
conditions  can  never  be  satisfied  exactly  with  a  finite  line  array  of 
discrete  transducers. 

Fortunately,  in  most  applications  It  Is  not  essential  that  the 
Kerr-Murdock  conditions  be  satisfied  perfectly,  as  will  be  demonstrated 
in  the  next  section.  Before  proceeding  we  observe  that,  due  to  the 
Integral  relationship  between  w^(x)  and  w^x),  their  Fourier  transforms 
are  related  by 

W,(f)  =  jZirfW^f)  ,  (30) 

and  their  moment  integrals  M1  n  and  n  obey 
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(31) 


The  latter  result  follows  easily  from  Integration  by  parts. 


VIII.  USING  A  GENERAL  ODD-EVEN  BEAMFORMER  PAIR 


We  will  consider  if  It  Is  possible  to  find  an  even  beamformer  other 
than  the  K-M  companion  beamformer  whose  output  is  proportional  to  the 
temporally  integrated  output  of  a  given  odd  beamformer  through  some 
angle  dependent,  real  valued,  scalar  factor  R.  This  proposition  can  be 
put  Into  equation  form,  (^(t^RqgU),  where  the  subscripts  1  and  2  denote 
the  odd  and  even  beamformers ,  respectively,  and  can  also  be  expressed  in 
terms  of  Fourier  transforms  as 

(j2trf)“1  Q-j  (f )  =  RQ2(f)  .  (32) 

With  the  aid  of  Eq.  (26),  this  equation  may  be  solved  for  R, 


The  use  of  Eq.  (30)  gives  the  equivalent  fom 


R 


sine  WKM^f') 


(34) 


where  f 'A(f/c)s1ne.  Clearly  the  only  way  that  R  can  depend  upon  the 
wave's  angle  but  not  Its  frequency  is  for  the  function  W2(f)  to  be  Iden¬ 
tical  to  W^(f)  or  a  scalar  multiple  of  it.  Since  equality  of  the  Fourier 
transforms  implies  equality  of  w2(x)  and  w^(x),  our  question  is  answered 
in  the  negative,  at  least  for  an  Incident  wave  of  unknown,  unrestricted 
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form.  No  beamformer,  nontrivial  ly  distinct  from  the  K-M  companion 
beamformer,  achieves  the  desired  result  for  all  Incident  waves. 

But  If  the  frequency  spectrum  of  the  Incident  wave  Is  restricted, 
then  Eq.  (34)  need  not  be  satisfied  at  all  values  of  f  and  e.  In  the 
Important  special  case  where  the  Incident  wave  Is  a  sinusoid  of  frequency 
fQ,  the  equality  Is  satisfied  automatically,  since  from  Eqs.  (33)  and 
(34) 


sine  WKM 


WKM(f)  Kl(0’V 
V2irr  '  J2irf0k2(e.i 


using  the  directivity  function  notation  defined  earlier.  The  scalar 
factor  R(e)  Is  clearly  an  odd-symmetric,  real  valued  function,  as  can  be 
seen  by  examining  each  of  the  functions  In  Eq.  (35). 

In  fact,  the  beamformer  outputs  will  be  related  by  a  scalar  factor 
even  In  the  nonslnusoldal  case,  provided  the  wave  is  incident  at  near¬ 
broadside  angle.  This  reasoning  follows  from  the  derivative  series 
representations  of  q-|(t)  and  q2(t), 


q2(t)  =  M2jQ  p(t)  +  ... 


(37) 


which,  after  Integrating  the  first,  combine  to  give  the  following  pro¬ 
portionality  law  that  Is  valid  for  small  e, 


The  angular  range  over  which  this  simple  proportionality  applies  can  be 
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determined  from  the  frequency  domain  requirement  of  Eq.  (34).  Let  A 
denote  the  maximum  argument,  for  which  WKM(f)/W2(f)  is  approximately 
constant,  l.e.,  if  |f|<A,  then  W)<M(f)/W2(f)  =  WKM(0)/W2(0)  -  Q/M2  Q* 
-M^  ^/M2  q  by  Eq.  (31).  Notice  that  the  size  of  A  depends  upon  how 
closely  the  even  beamformer  approximates  the  K-M  companion  beamformer, 
being  quite  large  If  the  approximation  is  good.  Now,  let  f  u  denote 
the  largest  frequency  for  which  the  Incident  wave  has  significant  spec¬ 
tral  energy.  Then  Eq.  (34)  admits  the  approximation 


over  the  angular  range  (sine |<cA/fmau.  Thus,  R  is  shown  to  be  angle 

nlaX 

dependent  but  frequency  independent  within  the  limits  of  this  approxima¬ 
tion. 


These  results  can  be  summarized  briefly.  A  proportionality  rule, 

q-|(t)  =  R(e)  q2(t)  ,  (40) 

holds  between  the  temporally  integrated  output  of  an  odd  beamformer  and 
the  output  of  an  even  beamformer,  subject  to  following  restrictions: 

Case  1 :  If  the  even  beamformer  is  the  K-M  companion  to  the  odd 
beamformer,  then  the  rule  holds  for  any  Incident  wave  and  for  all  e,  with 
R(e  )=ks=c_  1  sine. 

Case  2:  If  the  Incident  wave  is  sinusoidal  with  frequency  fQ,  the  rule 
holds  for  all  angles,  for  any  odd-even  pair  of  beamformers,  with  R(e) 
given  by  Eq.  (35). 

Case  3:  If  the  Incident  wave  is  nonsinusoidal ,  the  rule  still  holds  for 
any  odd-even  beamformer  pair,  but  only  for  angles  near  broadside. 


| e | <arcs1n(cA/ffmax) ,  where  the  parameters  A  and  fmax  are  defined  just 
prior  to  Eq.  (39), 

Once  R  Is  determined,  the  wave  angle  -Tr/2<e+ir/2  can  be 
unambiguously  determined  In  Case  1,  since  sine  Is  monotonic  over  that 
range.  The  same  principle  applies  to  Case  3  over  a  restricted  angular 
range  as  specified.  But  In  Case  2  the  dependency  of  R  upon  e  Is  com¬ 
plicated  by  the  nonconstant  ratio  term,  and  R(e)  will  probably  not  be 
monotonic  over  the  entire  range  -n/2<Q<+n/2.  Thus  the  unambiguous 
direction-finding  capability  will  be  limited  to  the  angular  range  over 
which  R(e)  is  monotonic.  However,  for  all  of  the  cases  listed  above,  It 
will  be  Impossible  to  determine  the  proportionality  factor  R  whenever 
the  wave  Is  incident  at  such  an  angle  that  the  even  beamformer  output  q2(t) 
Is  too  weak  to  be  measured  accurately.  Hence  the  angular  range  will  be 
constrained,  for  practical  purposes,  by  the  beam  pattern  of  the  even 
beamformer. 

Integrating  the  odd  beamformer  output  to  make  it  proportional  to 
the  even  beamformer  output  offers  a  possible  benefit  even  when  the  trans¬ 
form  method  Is  used.  The  transforms  Q-j (f )  and  Q2(f)  In  Eq.  (26)  will  be 
distorted  by  the  Inevitable  windowing  effects  of  a  computational  Fourier 
transform,  and  an  Imbalance  may  thus  result  due  to  the  dissimilarity  of 
q-j(t)  and  q2(t).  If  q-j(t)  is  Integrated  before  digital  sampling,  how¬ 
ever,  the  windowing  effects  will  be  balanced.  Thus  Q^f)  in  Eq.  (26) 
may  be  replaced  by  (j2trf (f) ,  where  (f)  denotes  the  computed 
Fourier  transform  of  q^t). 
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IX.  RELATED  PROCESSING  TECHNIQUES 

Those  familiar  with  radar  systems  will  have  noticed  that,  when 
restricted  to  sinusoidal  waves,  much  of  what  has  been  presented  here 
strongly  resembles  the  type  of  signal  processing  performed  by  a  monopulse 
radar.  The  same  processing  method  was  used  successfully  In  the  bearing 
deviation  Indicator  (BDI)  that  was  developed  at  Harvard  Underwater  Sound 

O 

Laboratory  during  World  War  II.  This  instrument,  which  sensed  the 
bearing  angle  of  sonar  echoes  from  a  hostile  submarine,  was  perfected 
prior  to  the  development  of  monopulse  radar.  Although  long  since 
declassified,  the  BDI  project  received  little  attention  In  the  open 
literature.  In  sharp  contrast,  monopulse  radar  has  been  elaborated  and 
extended  by  numerous  papers9  and  at  least  one  book.^ 

Among  other  things,  our  treatment  extends  to  general,  nonsinusoidal 
waves.  This  may  not  seem  pertinent  to  many  applications,  for  example  to 
active  sonar  systems  that  employ  sinusoidal  "pings,"  however,  It  is 
actually  Important  in  even  that  case  when  one  must  ponder  the  effects 
upon  bandwidth  of  using  short  pulse  durations  and  tuned  filters.  It  is 
significant  that  an  integrator,  as  we  have  shown,  provides  the  proper 
mechanism  to  extend  the  traditional  90°  phase  shift  operation  to  broadband 
signals  (rather  than  a  Hilbert  transform,  for  example),  when  coupled  with 
the  use  of  odd-even  shading  functions  that  satisfy  the  Kerr-Murdock 
conditions.  To  the  author's  knowledge,  the  radar  literature  contains  no 
references  to  the  use  of  the  Kerr-Murdock  conditions  for  anything  other 
than  providing  a  "linear"  (with  respect  to  sine)  sensed  angle  output,  for 
waves  of  a  single  frequency. 

Although  restricted  to  sinusoidal  waves,  the  monopulse  literature 
Includes  error  effect  studies,10  statistical  analyses,11  performance  data,12 
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and  design  criteria  that  may  be  of  considerable  value  in  predicting  the 
performance  of  paired  beamformers  used  in  direction  finding.  To  acquaint 
the  reader  with  monopulse  techniques  and  terminology,  and  to  make  the 
monopulse  literature  more  accessible,  the  following  detailed  summary  is 
provided.  The  reader  Is  hereby  cautioned  that  transitions  between  the 
sonar  and  radar  environments  will  be  made  without  comment. 

The  term  monopul se^  is  used  to  describe  radar  and  radio  tracking 
systems  that  measure  a  target's  azimuth  and/or  elevation  angle  after 
receiving  just  one  pulse  of  radio  frequency  (RF)  energy,  in  contrast  to 
systems  that  scan  the  target  by  slowly  rotating  the  antenna,  by  wobbling 
some  component  of  the  antenna,  or  by  sequentially  switching  between 
elements  of  a  compound  antenna.  In  the  monopulse  method  separate  beams 
are  used.  They  are  either  pointed  in  the  same  direction  with  displaced 
phase  centers  ("phase  comparison"  or  "Interferometric") ,  or  aimed  in 
slightly  different  directions  separated  by  a  squint  angle  ("amplitude 
comparison") .  The  dual  beams  may  be  obtained  from  paired  antennas,  or 
by  merely  employing  two  RF  feedpoints. 

Monopulse  radar  was  developed  primarily  to  circumvent  a  fault  of 
scanners,  namely  that  modulations  produced  by  scanning  are  similar  to  the 
modulations  produced  by  target  fluctuations,  resulting  in  possible 
confusion  or  Interference.  But  the  monopulse  method  quickly  gained  a 
reputation  for  remarkably  precise  measurements  of  target  direction  under 
good  signal-to-noise  ratio  (S/N)  conditions.  The  technique  has  given 
single  target  tracking  performance  with  angular  precision  that  bettered 
the  Rayleigh  resolution  criterion  by  2000  to  1  (only  200  to  1  if  all 
system  errors  and  biases  were  Included.)  As  a  device  for  precision 
direction-finding,  the  monopulse  radar  is  the  most  primitive  member  of 
a  family  of  advanced  nonlinear  processing  methods  that  Includes  super- 
resolution  techniques,  product  arrays,  maximum  likelihood  and  maximum 
entropy  processing,  adaptive  null  steering,  and  autoregressive  or 
Prony-type  analysis  of  angular  spectra.  Yet  it  is  the  only  method  that 
has  found  such  widespread  use  in  operational  radar  systems  as  a  means  of 
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extending  angular  precision  beyond  the  classical  limits  of  linear  array 
processing. 

Since  digital  computers  were  not  yet  available  when  the  monopulse 
method  was  developed,  estimates  of  target  angle  had  to  be  computed  with 
simple  electronic  circuitry,  and  that  fact  largely  determined  which 
statistic  would  be  computed.  Most  monopulse  systems  were  designed  to 
compute  the  real  part  of  the  ratio  of  the  complex  envelopes  of  the 
difference  and  sum  of  the  outputs  of  a  pair  of  slightly  displaced  beams 
(the  difference  signal  having  been  delayed  by  90°  In  the  case  of  phase 
comparison  monopulse).  This  quantity  is  proportional  to  the  target  angle 
over  a  limited  range. 

Although  monopulse  radars  typically  measure  both  azimuth  and 
elevation  angles,  for  simplicity  let  us  consider  only  the  measurement  of 
a  single  target  angle  e,  and  assume  that  the  target  is  confined  to  the 
same  plane  that  contains  a  pair  of  identical  line  arrays  used  to  observe 
It.  For  phase  comparison  monopulse,  these  receiving  subarrays  are 
col 11  near  but  their  centers  are  displaced  by  a  small  distance  2h.  For 
amplitude  comparison  monopulse,  they  have  a  common  center  but  their  axes 
are  separated  by  a  small  angle  28,.,  where  e  is  called  the  squint  semiangle 

w  S> 

(l.e.,  half  the  squint  angle  that  separates  the  beams).  The  left  and 
right  outputs  q^(t)  and  qR( t )  are  formed  by  identical,  real -symmetric 
shading  functions.  The  polar  coordinate  axis  for  measuring  target 
angle  0  is  oriented  as  shown  in  Fig.  1.  As  a  consequence  of  this  orienta¬ 
tion,  qL(t)=qR(t)  whenever  the  target  is  located  at  0=0. 

Recall  from  Eq.  (16)  that  the  directivity  function  of  a  line  array 
beamformer,  K(e,f0),  expresses  the  amplitude  and  phase  relationship  between 
the  sinusoidal  field  at  the  center  of  an  array  p(t)  and  its  output  q(t), 

Q  -  K(e,f0)P  ,  (41) 
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In  which  P  and  Q  arc  the  phasors  associated  with  p(t)  and  q(t),  and  6 
Is  the  target  angle  with  respect  to  the  array's  broadside  aspect.  If 
this  equation  were  applied  without  modification  to  both  component 
subarrays  of  a  monopulse  system,  an  Incompatibility  In  the  coordinate 
reference  would  result.  In  phase  comparison  monopulse,  the  subarray 
centers  are  not  coincident,  and  therefore  the  P's  would  be  different.  In 
amplitude  comparison  monopulse,  the  centers  are  coincident,  but  the 
coordinate  axes  for  d  would  disagree  due  to  the  angular  separation.  The 
necessary  modification  Is  simple.  Let  K(e)  denote  the  directivity 
function  that  both  subarrays  would  have  if  they  were  hypothetically 
coincident  (l.e..  If  h=0  or  eh=*0).  We  have  suppressed  the  notatlonal 
dependence  of  K  upon  f  for  convenience.  Recall  from  Eq.  (12)  that 
K ( e ) *=W[ ( f 0/c ) s  1  n e] ,  where  W( f )  Is  the  real,  even-symmetric  Fourier  trans¬ 
form  of  w(x),  the  aperture  shading  function  common  to  the  hypothetically 
coincident  subarrays.  But  suppose  we  continue  to  let  P  represent  the 
phasor  associated  with  the  field  p(t)  at  the  coordinate  origin,  even 
after  the  subarrays  are  separated  or  squinted.  Then  if  QL  and  QR  are 
the  phasors  associated  with  the  subarray  outputs  qL(t)  and  qR(t),  we  have 
for  amplitude  comparison  monopulse 

Ql  3  K(e-eh)P  =  K(eh-e)P  ,  (42) 

Qr  =  K(eh+e)P  ,  (43) 

while  for  phase  comparison  monopulse, 

Ql  =  ejB  K(e)P  ,  (44) 

QR  =  e"jB  K(e)P  ,  (45) 

where  B  k  2irf0c  h  sine,  which  Is  the  phase  shift  due  to  displacement  of 
the  subarray  center  by  distance  h. 

The  transformation  to  the  sum-and-difference  outputs 
qi (t)=q^(t)-qR(t)  and  q2(t)=qL(t)+qR(t)  is  obviously  Invertible,  so  these 
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outputs  can  be  used  for  all  subsequent  signal  processing  with  no  loss 
of  information.  If  Q-j  and  Q2  denote  their  phasors,  then  Q^*K^(e)P  and 
Q2"K2(0)P»  where 


K,(0)  4  Kl(0)  -  Kr(6) 

(46) 

K2(0)  4  KL(e)  +  KR(e) 

(47) 

For  phase  comparison  monopulse, 

K^e)  =  K(e)  ejB  -  K(e)  e"JB  =  2jK(e) 

sin8  , 

(48) 

K2(e)  »  K(e)  ejB  +  K(e)  e'jB  =  2K(e)  cosB 

(49) 

which  are  real-even  and  imaginary-odd,  respectively, 
unsurprising  since  the  displaced  subarrays  are  still 
system  may  be  regarded  as  a  single  line  array  with  a 
that  have  aperture  shading  functions 

These  results  are 

col  linear,  so  the 
pair  of  beamformers 

w-|  (x)  =  w(x+h)  -  w(x-h)  , 

(50) 

w2(x)  =  w(x+h)  +  w(x-h) 

(51) 

Since  w(x)  is  even,  w^x)  and  w2(x)  are  odd  and  even,  respectively.  One 
of  the  outputs  q-|(t)  or  q2(t)  must  be  sub<iected  t0  a  90°  phase  shift  to 
bring  them  into  phase  agreement.  If  we  use  the  symbol  J  to  represent 
the  90°  phase  shift  operator,  which  may  be  implemented  in  a  variety  of 
ways  for  sinusoidal  signals,  then  from  Eqs.  (48)  and  (49)  q-j (t)=JtanBq2(t) . 
Since  J2wf0q^(t)sq.(t)  for  a  sinusoid,  recalling  that  q^  denotes  the 
Integral  of  q^ ,  it  follows  that  a  phase  comparison  monopulse  system  with 
sum-and-difference  processing  constitutes  our  Case  2  odd-even  beamformer 
pair,  for  which  q^ (t)=R(e)q2(t).  Clearly, 
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R(e )  »  (2*^)“'*  tanB  *  (2^^)”^  tan(2irfohc^s1n0)  .  (52) 

For  this  case,  the  angle  can  always  be  determined  unambiguously  from  R(e) 
over  the  range  | s1ne|<c/4f0h. 

On  the  other  hand,  for  amplitude  comparison  monopulse 

K1 (e)  =  K(es-e)  -  K(es+e)  ,  (53) 

K2(e)  =  K(es-e)  +  K(es+e)  .  (54) 

These  directivity  functions  are  real -even  and  real -odd,  so  Instead  of 
being  90°  out  of  phase,  the  outputs  are  exactly  In  phase.  This  departure 
from  the  behavior  of  odd -even  beamformers  should  have  been  expected,  since 
the  component  subarrays  are  not  col  11  near.  However,  It  is  true  that 


where  the  quantity  In  brackets,  being  a  ratio  of  odd  and  even  real  functions, 
Is  Itself  an  odd  function,  monotonic  over  some  angular  range,  and  linear 
for  small  e.  Thus  e  can  be  determined  from  the  scalar  factor  that  relates 
and  q2,  without  the  need  for  Integrating. 

For  a  single  Incident  wave,  the  ratio  Q.|/Q2  (°r  the  case  of 

phase  comparison  radar)  ought  to  be  real.  However,  due  to  a  variety  of 
corruptive  effects,  It  may  possess  an  Imaginary  part.  It  is  customary  to 
use  the  real  part  of  the  ratio  for  angle  detection,  i.e. ,  to  use 


The  Implications  of  this  custom  will  be  discussed  In  a  later  section  of 
this  report. 
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Before  concluding  this  summary  of  monopulse  radar,  we  should  alert 
the  reader  to  the  fact  that  some  phase  comparison  systems  compare  the 
phases  of  q^(t)  and  q2(t)  directly  without  ever  forming  the  sum-and- 
dlfference  signals,  and  primitive  amplitude  comparison  systems  sometimes 
detect  the  amplitudes  separately  before  comparison.  Note  that  this  latter 
method  Is  not  strictly  equivalent  to  sum-and-dlfference  processing.  More¬ 
over,  the  literature  gives  considerable  attention  to  various  ways  of 
transforming  the  signals  prior  to  detection,  mostly  based  on  the  following 
two  simple  properties  of  complex  variables:  If  and  Z2  are  complex 
phasors  which  have  equal  amplitude  but  which  differ  In  phase  by  an  angle  e 
then  j(Z^+Z2)  and  Z^-Z2  have  like  phase  but  differ  In  amplitude  and  have 
the  ratio  tan(e/2).  On  the  other  hand,  If  Z-j  and  Z2  have  like  phase  but 
possess  an  amplitude  ratio  r,  then  +JZ2  and  Z-j-jZ2  have  like  amplitude 
but  differ  In  phase  by  2arctan(r). 
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X.  PHASE  SHIFT  SQUINTED,  AMPLITUDE  COMPARISON  MONOPULSE 
VERSUS  ODD-EVEN  BEAMFORMING 

In  a  sonar  implementation  of  an  amplitude  comparison  monopulse 
system,  It  Is  likely  that  the  beam  squinting  would  be  Implemented  by 
electronic  phase  shift  steering  of  a  single  line  array  rather  than  by 
using  two  physically  distinct  arrays.  For  this  case,  and  with  the  proper 
interpretation,  there  is  an  equivalence  to  the  use  of  an  odd-even  beam- 
former  pair. 

The  output  q(t)  of  a  simple  array  with  an  even-symmetric  aperture 
shading  function  w(x)  is,  from  Eq.  (1), 


sinej  dx 

If  the  complex  function  P(x)  expresses  the  amplitude/phase  of  the  inci¬ 
dent  sinusoidal  field  at  point  x,  and  Q  is  the  output  phasor,  then 


Q  =  f  w(x)  P(x)  dx  .  (57) 

The  array  can  be  steered  to  any  angle  es  by  applying  a  progressive 
phase  shift  along  the  array  sensor  axis: 

Q  -  J' w(x)  e“^Sx  P(x)  dx  , 

3  /  w(x)  c°s(Sx)  P(x)  dx  -  j  j' w(x)  sin(Sx)  P(x)  dx  ,  (58) 
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where  S»(f0/c)s1ne$.  If  an  amplitude  comparison  monopulse  system  Is  then 
formed  by  taking  sum-and-dlfference  outputs  and  Q2  after  phase  shift 
steering  to  the  two  squinting  angles  the  result  Is 


Ql  a  -if  w,(x)  P(x)  dx 

(59) 

Q 2  *  f  w2(x)  dx  » 

(60) 

where  the  odd-even  shading  functions  and  w2  are 

W'j  (x)  =  2w(x)  sin(|  fQ  slnej 

(61) 

w2(x)  =  2w(x)  cos(|  fQ  s1neh) 

(62) 

From  Eq.  (60)  we  see  that 

d2(t)  -/ w2(x)  p(t  -  ~  sine)  dx 

(63) 

A  similar  result  holds  for  Eq.  (59)  except  that  we  must  Insert  the  pre- 

A 

vlously  defined  symbol  J  to  account  for  the  90°  phase  shift  that  is 
Implied  by  the  j  factor: 


Jqi  (t)  -f  Wl(x)  p(t  -  |  sine)  dx 


(64) 


Examination  of  Eqs.  (63)  and  (64)  reveals  the  phase  shift  squinted, 
amplitude  comparison  monopulse  system  to  be  equivalent  to  an  odd-even  pair 
of  beamformers,  except  that  It  already  has  a  90°  phase  shift  In  Its  q-j ( t) 
output.  But  that  phase  shift  was  effected  artificially  as  part  of  the 
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phase  shift  steering  process,  and  could  just  as  well  have  been  postponed 
until  needed  to  bring  the  sum-and-dlfference  outputs  Into  phase  agree¬ 
ment.  Thus  the  difference  between  these  two  methods  Is  one  of  Inter¬ 
pretation.  The  situation  is  quite  different  with  a  true  amplitude 
comparison  system,  which  performs  its  90°  phase  shift  by  spatial  pro¬ 
cessing  rather  than  electronic  or  computational  processing. 
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XI.  THE  PARTICLE  MOTION  EQUIVALENCE  PRINCIPLE 


The  output  q2(t)  of  an  even  beamformer  Is  related  to  the  incident 
wave  by  Eq.  (4),  q2(t)=w2  0(t)*p(t),  where  w2>e(t)=  ^w2(t/k),  with 
k=sine/c.  Suppose  the  even-symmetric  shading  fuction  w2(x)  has  been 
selected  to  yield  a  highly  directive  broadside  beam.  Then  W2(f)  must 
mimic  the  frequency  response  of  a  low  pass  filter  if  the  directivity 
function  K2(e,f)=W2[(f/c)sine3  is  to  pass  only  those  waves  having  small 
0.  Equation  (6)  gives  the  frequency  domain  behavior, 
Q2(f)=W2[(f/c)sine]P(f),  by  which  the  beamformer  may  be  characterized  as 
an  angle  dependent  low  pass  filter  whose  cutoff  frequency  is  inversely 
proportional  to  sine,  or  as  an  angular  filter  whose  small  angle  cutoff, 
i.e.,  beamwidth,  is  inversely  proportional  to  frequency.  This  spatial 
filter  Is  remarkable  in  being  devoid  of  phase  shift  at  all  frequencies, 
since  W(f)  takes  only  real  values.  Let  fL  denote  the  low  pass  cutoff 
frequency,  so  that  |W(f)-W(0)[<<5  for  |f|<f^,  where  the  tolerance  parameter 
6  may  be  chosen  to  fit  the  3  dB  point,  for  example.  We  define  broadside 
window  as  the  set  of  frequency  and  angle  combinations  for  which 

£  sine  <  f^  .  (66) 

Thus,  if  an  incident  wave  "passes  through"  this  window,  i.e.,  if  all  of 
its  spectral  energy  is  at  frequencies  that  satisfy  Eq.  (66),  then 
Q2(f)  fa  W2[(f/c)sine]P(f)  =  W2(0)  P(f)  and  it  follows  that 

q2(t)  =  W2(0)  p(t)  .  (67) 

(For  simplicity  we  shall  drop  the  approximation  symbol.)  Figure  4 
illustrates  the  windowing  effect  of  the  even  beamformer's  directivity 
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SHIELDING  EFFECT  OF  THE  EVEN  BEAMFORMER 
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function  upon  waves  coming  from  sources  at  various  angles.  The  effect 
may  be  regarded  as  a  shielding  action.  This  figure  must  be  Interpreted 
with  caution  since  the  opacities  of  the  directional  response  are 
represented  only  figuratively,  with  no  frequency  dependence  shown. 

Now  suppose  this  even  beamformer  Is  combined  with  its  K-M  companion 
beamformer  to  form  an  odd-even  pair.  Also,  for  the  moment,  let  us 
restrict  our  discussion  to  "plane"  waves  emanating  from  very  distant 
wave  centers  that  lie  in  a  single  X-Y  plane  that  contains  the  receiving 
aperture.  If  a ( t )  denotes  the  acoustic  particle  displacement  along  the 
ray  path  of  a  solitary  incident  plane  wave,  then  p(t)=-(B/c)a(t) ,  where 
B  denotes  the  bulk  modulus.  The  scaling  of  the  aperture  shading  function 
w2(x)  Is  arbitrary,  so  assume  for  mathematical  convenience  that  it  Is 
scaled  to  make  W2(0)=-(c/B) .  With  this  calibration  adjustment  the  beam- 
former  output  for  a  solitary  plane  wave  incident  at  angle  e  is  given  by 
q2(t)aa(t),  provided  the  wave  passes  through  the  broadside  window 
unscathed.  (Note  that  we  are  assuming  a  line  array  of  acoustic  pressure 
sensors  that  accurately  measure  the  Instantaneous  pressure  of  any  plane 
wave  arriving  within  its  broadside  window.)  It  follows  that  the  particle 
displacement  could  be  determined  merely  by  Integrating  q2(t),  i.e., 

a ( t )  =  q2(t)  .  (68) 

To  represent  the  acoustic  particle  motion  on  an  X-Y  display  device 
such  as  an  oscilloscope,  it  would  be  necessary  to  set  Y=a(t)cose  and 
X*-a(t)sin0.  Recall  that 

q,(t)  *  2Mq2(t)  . 

so  we  may  simply  set 

X  =  -cq1 ( t )  ,  (69) 
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where  the  double  bar  denotes  two  temporal  integrations.  Unfortunately, 
there  Is  no  linear  way  to  form  Y=a(t)ccse  without  prior  knowledge  of  e, 
but  if  the  broadside  window  admits  only  waves  whose  angles  are  small 
enough  that  cose  is  approximately  one,  then  we  can  set 

Y  =  q2(t)  .  (70) 

Notice  that  the  error  in  setting  cose=l  is  at  most  10X  for  angles  less 
than  25°,  or  IX  for  angles  less  than  8°,  so  the  distortions  in  the  dis¬ 
play  will  be  quite  small  in  most  cases  of  interest. 

For  a  solitary  incident  wave,  the  particle  motion  (actual  or 
simulated)  occurs  along  a  straight  line  that  points  toward  the  wave  cen¬ 
ter,  presumed  to  be  an  infinite  distance  away.  But  this  display  con¬ 
tinues  to  accurately  represent  the  particle  motion  at  the  point  of 
measurement  even  when  several  distinct  waves  with  arbitrary  spectra 
arrive  from  a  variety  of  angles  near  broadside.  This  is  true  because 
the  process  of  beamforming,  scaling,  and  integration  is  just  as  subject 
to  linear  superposition  as  is  the  acoustic  medium  itself.  Note  also 
that  the  integrations  are  reversible,  so  the  display  preserves  all  of 
the  information  contained  in  the  beamformer  outputs.  The  eye  will  likely 
not  be  able  to  interpret  that  information  as  well  as  an  electronic  angle 
detector  of  good  design.  The  subject  of  angle  or  slope  detection  will 
be  covered  later.  For  the  moment  we  simply  note  that,  loosely  speaking, 
the  angle  of  a  dominant  wave  in  a  background  of  weaker  waves  is  revealed 
by  the  predominant  direction  of  particle  motion. 

Any  odd-even  beamformer  pair  may  also  be  used  to  generate  the 
display,  provided  the  incident  waves  have  acceptable  frequency  spectra 
and  angles  sufficiently  near  broadside  so  that  q-j  (t)=R(e)q2(t),  and  so 
that  R(e)  accurately  approximates  (sine)/c  to  within  a  multiplicative 
factor.  The  X  variable  must  be  properly  scaled  to  compensate  for  that 
factor,  of  course.  Moreover,  it  may  be  possible  to  design  a  substitute 
beamformer  whose  output  is  related  to  that  of  the  original  even 
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beamformer  by  the  factor  cose  over  a  limited  frequency  range.  This  would 
reduce  the  distortion  that  results  from  using  q2(t)  to  drive  the  Y 
variable. 

The  significance  of  the  argument  just  presented  goes  beyond  the 
mere  generation  of  an  interesting  display.  It  demonstrates  that  the 
beamformer  outputs  contain  no  more  and  no  less  information  than  is 
required  to  describe  the  acoustic  particle  motion  at  a  hypothetical 
measurement  point  that  is  shielded  by  the  directivity  function  associated 
with  the  even  beamformer.  Nothing  can  be  determined  about  the  individual 
waves,  or  about  the  locations  of  their  wave  centers,  that  does  not  reveal 
itself  in  the  acoustic  motion  of  this  single  particle.  The  problem  of 
properly  interpreting  the  beamformer  outputs  can  be  approached  purely 
from  the  standpoint  of  basic  physics  and  signal  geometry,  with  little 
regard  to  the  parameters  of  the  line  array  that  might  be  used  to  obtain 
them. 


The  particle  motion  equivalence  principle  just  stated  clarifies  the 
role  of  directivity  in  identifying  the  geometrical  configuration  of  a 
distant  cluster  of  wave  centers.  The  direction  of  a  solitary  wave  cen¬ 
ter  can  always  be  determined  with  infinite  precision  in  the  absence  of 
measurement  noise.  But  if  waves  with  similar  frequency  spectra  arrive 
from  a  cluster  of  wave  centers  whose  angular  extent  falls  well  within  the 

broadside  window,  then  it  may  not  be  possible  to  separate  their  effects 

upon  particle  motion.  Nevertheless,  analysis  of  its  motion  may  provide 

some  information  about  the  cluster,  and  it  would  seem  that  such  knowledge 

could  only  be  enhanced  by  sharpening  the  directivity  function.  But 
unless  the  directivity  can  be  improved  to  such  a  degree  that  the  broad¬ 
side  window  masks  some  of  the  elements  of  the  cluster,  there  will  be  no 
effect  on  the  particle  motion,  and  hence  no  effect  on  the  information 
content  of  the  beamformer  outputs.  This  seems  to  suggest  that  a  modest 
enlargement  of  the  receiving  aperture,  although  it  would  certainly 
improve  the  directivity,  might  be  of  no  aid  in  identifying  the  cluster 
structure.  Yet,  under  some  circumstances,  that  is  definitely  not  true. 

If  the  beamformer  outputs  are  seriously  corrupted  by  measurement  noise 


47 


of  a  sort  that  does  not  grow  with  the  size  of  the  receiving  aperture, 
such  as  electrical  noise,  then  Its  enlargement  should  increase  the 
received  signal  energy,  thereby  Improving  S/N  even  though  the  enlarge¬ 
ment  has  no  effect  on  the  particle  motion.  This  improvement  would  cer¬ 
tainly  be  beneficial  for  determining  the  direction  of  a  solitary  wave 
center,  since  In  that  instance  the  precision  Is  limited  only  by  measure¬ 
ment  noise.  But  in  the  case  of  waves  from  a  cluster,  the  analysis  of 
particle  motion  might  be  limited  more  by  Its  sheer  complexity  than  by 
measurement  noise,  and  a  modest  aperture  enlargement  would  then  be  of  no 
benefit. 

It  should  be  acknowledged  that  this  presentation  has  ignored  some 
acoustical  details.  We  have  assumed  perfect  transduction,  and  the  dis¬ 
tortion  of  particle  motion  caused  by  the  physical  presence  of  the  array 
has  not  been  considered.  The  particle  motion  discussed  here  is  predicated 
upon  hypothetical  free  field  conditions,  i.e.,  it  is  the  motion  at  the 
measurement  point  in  the  absence  of  the  array,  but  shielded  by  the  even 
beamformer's  directivity  function.  We  have  also  used  the  terms  plane 
wave  and  wave  center  despite  their  contradictory  nature.  This  practice 
will  be  continued,  under  the  assumption  that  the  wave  centers  are  so  dis¬ 
tant  that  the  incident  waves  are  effectively  planar.  Notice  that  setting 

X=-cq1(t)  (71) 

and 

Y  =  q2(t)  (72) 

provides  a  display  of  particle  velocity,  which  is  just  as  useful  for 
determining  wave  direction.  For  this  reason,  the  beamformer  outputs 
qij(t)  and  q2(t)  may  be  referred  to  as  the  velocity  outputs. 

The  use  of  an  oscilloscope  display  to  reveal  wave  angle  (using  sum- 

O 

and-difference  beamformers)  has  been  known  since  the  1940's,  and  has  been 
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discussed  In  at  least  one  textbook.  But  In  Its  traditional  form  the 
display  has  a  frequency  dependent  angular  scale,  and  the  connection 
between  the  display  and  acoustic  particle  motion  has  been  Ignored. 


XII.  ANALYSIS  OF  X-Y  DISPLAY  MOTION 
FOR  MULTIPLE  WAVES 


Regardless  of  whether  one  constructs  a  particle  displacement  display 
or  a  particle  velocity  display,  or  merely  connects  the  q^(t)  and  qg( <:) 
outputs  of  any  odd-even  beamformer  pair  to  an  X-Y  display  device  with 
arbitrary  scaling,  a  solitary  Incident  wave  will  manifest  Itself  as  a 
straight  line  through  the  origin.  If  the  display  point  (X,Y)  Is  (q^qg), 
then  q‘-j{t)=R(e)q2(t),  with  R(© )  given  by  Eq.  (35),  and  the  point  lies  on 
the  line  X-RY«0  at  all  times.  The  displayed  point  will  trace  out  that 
line  repeatedly  with  each  oscillation  of  the  Incident  wave.  The  line 
will  tilt  to  the  right  or  left  of  vertical  in  accordance  with  the  mapping 
of  Incidence  angle  to  display  angle  that  Is  induced  by  the  function  R(e) 
(see  Fig.  5).  The  length  of  the  line  as  projected  onto  the  Y-axis 
represents  the  peak-to-peak  amplitude  of  the  output  of  the  even  beam- 
former,  so  It  will  be  relatively  Independent  of  the  angle  of  incidence 
so  long  as  the  ray  path  remains  within  the  broadside  window. 

If  two  waves  arrive  at  different  times  and  different  angles,  then 
two  lines  will  be  traced  In  succession  (see  Fig.  6(a)).  If  (X^.Y^) 
denotes  the  coordinates  of  the  displayed  point  during  the  arrival  of  the 
first  wave  and  (X2,Y2)  denotes  those  during  the  second,  then  X^-RjY^O 
and  X2-R2Y2=0  where  R^  and  R2  denote  R(e-|)  and  R(e2).  But  If  the  waves 
arrive  simultaneously,  yielding  the  superposition  of  coordinates 
(XS,YS)»(X1+X2,Y1+Y2),  then  the  display  point  will  not  be  confined  to 
a  line.  Its  coordinates  (XS,YS)  will  obey 

Xs  -  R,YS  =  (i-R,R;’)x2  (73) 

and 
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FIGURE  6 

BOUNDARY  PARALLELOGRAM 


(74) 


xs  -  R2ys  .  (l-RjOx,  • 

and  if  we  assume  the  subscripts  have  been  chosen  so  that 

(l -R1  Rg1  )min{X2)  s  X$  -  R1V$  S  (l-R^j^maxiXg} 

and 

(RJ,R-1-l)m1n{-Xi }  s  Xs  -  R2YX  s  (r^’-iJmxI-X,) 

These  limits  constrain  the  display  point  (Xg,Ys)  to  the  interior  of  a 
parallelogram.  Two  of  its  sides  are  parallel  to  the  line  X-R1Y=0  with 

V  intercepts  Y^l-R^^nlX.,}  end  Y+*(l- R^mexlX.,),  and  the  other  two 
sides  are  parallel  to  the  line  X-R2Y=0  with  Y  intercepts 

Y  «(R2R^1-l)min{-X1)  and  Y+=(R2R^1-l)max{-X1 }.  If  wave  1  were  extinguished 
then  both  bounds  of  Eq.  (76)  would  become  equal ,  and  the  display  would 
collapse  to  a  single  line  segment  through  the  origin  with  slope  R2,  its 
ends  touching  the  parallogram  boundaries.  Similarly,  if  wave  2  were 
extinguished,  the  result  would  be  a  line  with  slope  R-j ,  also  touching 

the  parallelogram  boundaries  (see  Fig.  6(b)).  These  line  segments  bisect 
the  sides  of  the  parallelogram,  and  may  be  regarded  as  non-orthogonal  but 
additive  components  of  the  display  motion  (see  Fig.  6(c)). 

If  wave  2  is  periodic  with  period  T2,  then  X2  attains  its  maximum 
and  minimum  values  at  least  once  during  each  period,  i.e.,  the  display 
point  touches  both  sides  of  the  parallelogram  defined  by  Eq.  (75)  during 
every  time  interval  of  duration  T2-  Similarly,  if  wave  1  has  period  Tr 
then  the  remaining  two  sides  are  touched  at  least  once  during  that 
period.  If  both  waves  are  periodic,  then  unless  there  is  some  coherence 
between  the  two  waves,  it  is  likely  that  enough  different  points  of  the 
bounding  parallelogram  will  be  touched  to  define  it  unambiguously  after 
a  brief  interval  of  observation.  Then  R-j  and  R2  can  be  determined  from 
the  slopes  of  its  sides,  and  if  the  angles  of  incidence  lie  within  the 
range  where  R(e)  is  monotonic,  then  9-|  and  e2  can  be  determined  from  R-j 


then 
.  (75) 


(76) 
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and  R2  (see  Fig.  7).  It  Is  thus  possible  to  resolve  both  angles  of 
Incidence  when  the  Incident  waves  are  periodic  but  not  coherent,  even 
when  they  arrive  simultaneously.  But  It  might  be  preferable  just  to 
separate  them  by  using  tuned  filters. 

More  Important  Is  the  case  In  which  the  two  waves  are  coherent, 
being  sinusoids  of  the  same  frequency.  In  this  case,  tuned  filters 
would  be  of  no  use  In  separating  the  components.  Both  coordinates  X$(t) 
and  Y,j(t)  are  sums  of  sinusoids  of  the  same  frequency;  hence  both  reduce 
to  single  sinusoids,  not  necessarily  of  the  same  phase  or  amplitude. 
Anyone  familiar  with  Llssajous  patterns  knows  that  the  resulting  display 
will  be  an  ellipse.  As  a  result  of  the  properties  discussed  In  the  pre¬ 
ceding  paragraph,  the  ellipse  must  touch  all  four  of  the  parallelogram's 
sides,  l.e.,  It  Is  Inscribed  in  the  parallelogram  (see  Fig.  8(a)).  In 
the  special  case  where  the  waves  arrive  In  phase  the  ellipse  degenerates 
to  a  straight  line,  which  must  be  one  of  the  diagonals  of  the  parallelo¬ 
gram  (in  order  to  touch  all  four  sides).  The  long  and  short  diagonals 
thus  represent  constructive  and  destructive  Interference  of  the  waves. 
Unfortunately,  if  the  phase  relationship  of  the  incident  waves  remains 
steady,  the  ellipse  will  remain  steaJy  as  well.  A  single  inscribed 
ellipse  is  not  sufficient  to  define  the  parallelogram.  But  if  the  rela¬ 
tive  phase  between  the  wave  components  drifts  while  their  amplitudes 
remain  constant,  then  two  or  more  inscribed  ellipses  will  be  seen  and  the 
boundary  parallelogram  will  be  revealed,  i.e.,  the  two  wave  angles  will 
be  resolved  (see  Fig.  8(b),  (c),  (d)).  If,  as  in  active  sonar,  the  two 
waves  are  actually  echoes  from  objects  that  are  physically  separate, 
then  their  relative  phase  may  be  altered  by  physical  motion  of  the  trans¬ 
mitter,  the  receiver,  or  the  objects  themselves,  or  by  varying  the 
operating  frequency. 

The  inscribed  ellipse  interpretation  has  important  implications  for 
detecting  the  angle  of  a  supposedly  solitary  wave  that  has  been  cor¬ 
rupted  in  some  way.  Suppose  that,  due  to  multipath  (refractive  or 
reflective)  effects,  what  started  out  as  a  single  wave  arrives  as  two 
waves,  at  slightly  different  angles.  Or,  suppose  that  a  weaker  wave  of 
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FIGURE  7 

RESOLUTION  OF  DIRECTIONAL  COMPONENTS 
FOR  WAVES  OF  DIFFERENT  FREQUENCIES 
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the  same  frequency  arrives  from  a  different  wave  center,  perhaps  due  to 
reflection  of  a  transmitted  wave  from  a  boundary.  Given  that  the 
stronger  wave  produces  a  longer  line  segment  component  of  the  display 
than  the  weaker  wave,  it  is  clear  (see  Fig.  9)  that  the  apparent  direc¬ 
tion  of  the  wave,  as  interpreted  by  any  reasonable  characterization  of 
the  ellipse's  angle  of  tilt,  always  lies  within  the  longer  wedge  formed 
by  the  parallelogram  diagonals.  Thus  the  effect  of  the  weaker  wave  is 
either  to  attract  the  apparent  direction  towards  itself  if  the  inter¬ 
ference  is  constructive,  or  to  repel  it  if  the  interference  is  destruc¬ 
tive.  In  the  latter  case,  the  apparent  direction  does  not  lie  between 
the  actual  wave  centers.  This  phenomenon,  whereby  two  wave  centers  con¬ 
spire  to  appear  elsewhere,  is  known  as  "glint"  in  the  radar  monopulse 
literature.  Nevertheless,  if  the  offending  wave  is  significantly  weaker, 
for  example,  if  its  contributed  line  segment  component  of  the  display 
is  less  than  one  third  the  length  of  that  of  the  main  wave,  then  the 
angular  extent  of  the  "wedge  of  uncertainty"  is  relatively  small.  In 
any  case,  parallelogram  sketches  are  useful  tools  for  predicting  such 
effects.  An  even  clearer  picture  of  performance  can  be  obtained  from  an 
X-Y  oscilloscope  with  two  oscillators  set  at  slightly  differing  fre¬ 
quencies  to  simulate  waves  having  a  slowly  varying  phase  difference.  A 
means  must  be  provided  to  additively  mix  the  two  signals  into  both  chan¬ 
nels  of  the  oscilloscope. 

If  more  than  two  waves  arrive  at  the  same  time,  the  situation  becomes 
rather  more  complicated.  It  is  straightforward  to  show  that  the  gen¬ 
eralization  of  Eqs.  (75)  and  (76)  to  N  waves  is 

An  s  xs  -  Vs  s  Bn  for  n=1’2’-*-’N  »  (77) 

where  A  and  B  are  the  minimum  and  maximum  values  of 
n  n 

ntfn 
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EFFECT  OF  A  WEAKER  COMPONENT 


ARL.'UT 

AS-82-789 

TLH-GA 

6-21-82 


:  tf.  Witt  £Ah  Ah  t  k 


These  Inequalities  confine  the  display  point  to  the  Interior  of  a  polygon 
having  2N  sides,  grouped  as  N  pairs  of  parallel  sides.  Unfortunately 
for  N>2,  the  quantities  An  and  Bn  depend  upon  more  than  one  wave,  and 
even  if  each  wave  Is  periodic  there  is  no  guarantee  that  the  display 
point  will  touch  the  corresponding  sides  during  every  period.  But  for 
every  pair  of  sides  of  the  boundary  polygon  that  can  be  discerned,  the 
corresponding  slope  Rn  can  be  determined  and  the  direction  of  that  wave 
component  resolved.  The  feasibility  of  such  an  approach  for  a  given  set 
of  multiple  waves  is  best  explored  via  the  oscilloscope  method,  or  by  a 
computer  simulation. 

A  related  method  for  resolving  two  monopulse  radar  targets  has  been 
developed  by  Sherman, ^  but  that  method  is  purely  algebraic,  and  lacks 
the  benefits  of  a  geometrical  interpretation. 

Several  investigators  have  analyzed  monopulse  radar  behavior  in  the 
presence  of  multiple  targets,  employing  statistical  methods  in  most 
cases.  Recall  from  Eq.  (56)  that  a  monopulse  radar  computes 


r  =  Re 


or 


Re 


JQ 


2j 


If  the  output  phasors  and  Qg  represent  the  combined  radar  returns  from 
many  point  scatterers,  regarded  as  ideal  point  targets, 


«2* 


?  Q2.1 


(78) 


then 


r 


Re 


e  Q 


2,1 


(79) 
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where 


R 


i 


or 


Note  that  each  point  target  component  would  produce  a  real -valued  ratio 
that  accurately  indicates  its  angle  relative  to  broadside,  i.e., 
R-jaR(8j).  So, 


r 


r  R(0,)Q2>l 


r  Q 


2,1 


(80) 


Equation  (80)  may  be  used  as  the  starting  point  for  investigating  the 
behavior  of  monopulse-type  receivers  for  multiple  targets,  either  via 
statistical  methods  or  simulation.^0 
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XIII. 


SLOPE  DETECTION 


Regardless  of  whether  one  really  wishes  to  generate  it,  the  concept 
of  the  X-Y  display  suggests  good  ways  to  extract  wave  angle  from  the  odd- 
even  beamformer  outputs.  As  shown  in  the  previous  sections  of  this 
report,  the  beamformer  outputs  can  be  used  to  generate  signals  X(t)  and 
Y(t)  that,  for  a  solitary  Incident  wave  in  a  perfect,  noise-free  environ¬ 
ment,  obey  the  straight  line  relationship  X(t)=R(e)Y(t)  at  all  times,  and 
e  can  be  determined  from  R  within  some  monotonic  range.  We  may  even 
regard  the  observed  factor  R,  the  slope  of  the  displayed  line,  to  be  a 
direct  measure  of  wave  angle  in  nonlinear  "R  units".  (Although  R  will 
be  referred  to  as  the  "slope",  note  that  it  is  the  slope  of  X  versus  Y 
rather  than  the  more  familiar  Y  versus  X.)  Noise,  multipath,  and 
extraneous  wave  sources  may  distort  the  line  into  something  much  more 
complex,  such  as  an  ellipse  or  a  Brownian  motion  path.  Neverthless,  it 
is  to  be  expected  that  the  straight  line  component  of  the  dominant  wave 
will  assert  itself  as  a  line  of  elongation,  an  axis  of  symmetry,  or  some 
other  attribute  that  gives  the  display  an  identifiable  slope.  The  prob¬ 
lem  of  slope  detection  may  be  posed  in  a  gross  way  as  finding  an  estimate 

A 

R,  such  that 


X(t)  =  RY(t)  .  (81) 

If  9=0,  then  under  ideal  conditions  the  display  must  be  a  vertical  line 
with  R=0.  However,  in  the  presence  of  random  disturbances,  the  display 
will  distribute  itself  symmetrically  about  the  vertical  axis,  at  least 
in  a  statistical  sense,  so  that  the  product  X(t)Y(t)  averages  to  zero. 

If,  on  the  other  hand,  the  wave  angle  deviates  from  the  broadside  direc¬ 
tion  or,  equivalently,  if  the  line  array  is  physically  rotated,  then  the 
idealized  straight  line  will  have  a  nonzero  slope  R,  i.e.,  it  will  rotate 
through  an  angle  0,  where  R=tane.  If  the  wave  source  is  not  so  simple 
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and  has  an  angular  extent  that  produces  a  more  complex  display  than  a 
mere  line,  then  the  entire  display  pattern  will  rotate.  There  will  be 
little  or  no  discernible  distortion  associated  with  this  rotation  If  the 
odd  symmetric  function  R(e),  besides  being  monotonic,  Is  approximately 
linear  over  an  angular  span  that  encompasses  all  of  the  wave  sources 
which  contribute  significantly  to  the  display  pattern.  Moreover,  If 
statistically  independent,  Gaussian  background  noise  corrupts  the  display 
variables  X(t)  and  Y(t),  and  if  these  variables  have  been  scaled  for 
equal  noise  variance,  then  the  display  components  contributed  by  this 
noise  are  statistically  invariant  with  respect  to  rotation.  Thus,  under 
a  broad  class  of  circumstances,  we  are  led  to  estimate  the  slope  as 


A 

R  =  tan$  ,  (82) 

where  3  is  the  clockwise  angle  through  which  the  display  pattern  must  be 
rotated  to  decorrelate  the  coordinates  of  the  display  locus,  i.e.,  3  is 

V 

the  angle  required  to  make 


av{x( t)y (t ) }  ■  0  ,  (83) 


where 


x(t) 

y(t) 


"cos  3 
sin3 


sin3* 

cos3 


x 


(84) 


The  averaging  operator  denoted  by  av{*}  is  quite  general.  Our  subsequent 
developments  will  depend  upon  only  two  properties:  (1)  the  averaging 
process  must  be  a  linear  operation,  and  (2)  the  average  of  a  nonnegative 
quantity  must  not  be  negative.  Included,  for  example,  is  the  case  of  an 
"average"  formed  by  a  low  pass  filter. 


If  the  product  x(t)y(t)  is  formed  from  Eq.  (84)  and  then  substituted 
into  Eq.  (83),  and  the  result  divided  by  cos^3  and  av{Y^(t)>,  the  fol¬ 
lowing  equation  results: 


(85) 


2 

(1-tan  e)r  -  (tans)p  »  0 
where  r  and  p  are  defined  by 


r .  mimm  ,  m 

av{ Y^( t) > 


p  »  av{Y2(t)>  -  av(X2(t)> 
av{Y2(t)} 


(87) 


If  R  is  substituted  for  tang,  then  a  quadratic  equation  results: 


+ 


A 

pR 


r  =  0 


(88) 


whose  roots  are 


R, ,  R2  -  "P  1  ,1  iS  ,  (89) 

A 

where  R^  denotes  the  root  with  the  positively  signed  radical.  Note  that 
the  product  of  the  two  roots  is  -1,  i.e.,  each  is  the  negative  reciprocal 
of  the  other.  This  result  should  come  as  no  surprise  since  negatively 
reciprocal  slopes  imply  orthogonality,  confirming  that  the  rotated  dis¬ 
play  remains  just  as  uncorrelated  if  it  is  rotated  an  additional  90°. 

Of  these  two  orthogonal  choices  for  the  apparent  slope  of  the  display, 
it  is  logical  to  select  the  one  along  which  the  display  is  most  strongly 
distributed,  and  this  can  be  accomplished  by  picking  the  root  that  has 
the  same  sign  as  av{X ( t ) Y ( t ) } ,  i.e.,  the  same  sign  as  r.  But  the  radical 
term  in  Eq.  (89)  is  always  larger  in  magnitude  than  -p,  so  the  sign  of 
the  numerator  is  controlled  by  the  sign  in  front  of  the  radical.  Thus 
R-j  is  the  root  that  always  has  the  same  sign  as  r,  so  we  take 
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a 


2r 


(90) 


A  A 


R  ■  R,  “  -(R2) 


p  + 


The  quantity  In  brackets  must  lie  between  0  and  2*  as  can  be  shown  by 
some  mildly  tedious  algebra  that  will  only  be  outlined  here:  (1)  the 
radical  Is  larger  in  magnitude  than  p,  so  the  number  in  brackets  Is  non¬ 
negative.  (2)  From  the  two  properties  of  av{*}  given  above,  It  Is 

?  2  ? 
easily  proven  that  [av{XY>]  never  exceeds  the  product  av{X  }av{Y  },  by 

mimicking  the  standard  quadratic  discriminant  proof  from  basic  statistics, 

wherein  av{ • >  denotes  expectation.  (3)  It  follows,  after  some  simple 

2  2 

algebra,  that  the  expression  under  the  square  root  sign,  p  +4r  ,  is 
bounded  by  [(av{x2}+av{Y2})/av{Y2}]2.  (4)  The  sum  of  p  and  the  square 

root  of  this  bound  is  2,  so  the  quantity  in  brackets  does  not  exceed  1. 


A 

Hence,  the  estimated  slope  R  agrees  in  sign  with,  and  equals  or 
exceeds  in  magnitude,  the  coefficient  r.  Both  quantities  have  a  statis¬ 
tical  interpretation.  If  the  display  comprises  a  set  of  points  that 
would  fit  a  straight  line  through  the  origin  were  it  not  for  the  corrup¬ 
tive  effects  of  independent,  Gaussian  perturbations  of  zero  mean,  then  r 
is  the  sampled  coefficient  of  regression  of  X  upon  Y.  If,  in  addition, 
the  axes  have  been  scaled  for  equal  noise  variance,  then  R  is  the  maxi¬ 
mum  likelihood  estimate  of  slope.16 

When  the  display  is  dominated  by  waves  (including  directional 
noises)  arriving  at  angles  very  near  broadside,  the  X(t)  signal  will 
remain  very  weak,  since  it  is  derived  from  the  output  of  the  beamformer 
that  has  a  broadside  null.  Expressed  another  way,  the  display  will  be 
tall  and  thin.  Thus  r  will  be  quite  small  and  p  will  approximate  unity, 
so  that  r  is  an  approximate  estimate  of  R, 

r  2  R  =  R  (91) 


ryi~; 


If  the  beamformer  outputs  are  dominated  by  sinusoidal  waves  whose  spec¬ 
tral  energy  is  concentrated  at  some  frequency  fQ,  as  in  the  case  of 
active  sonar,  then  X ( t )  and  Y(t)  will  trace  an  ellipse,  degenerating  to 
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a  straight  line  for  a  solitary  Incident  wave  having  a  well  defined  angle. 

If  we  represent  the  amplitude  and  phase  of  the  sinusoids  X(t)  and  Y(t) 

by  phasors  Z  and  Z  ,  then  If  av{ • }  represents  a  long  term  arithmetic 
x  y 

mean  or  a  short  term  mean  over  an  Integral  number  of  periods,  then  we 
see  by  analogy  with  the  well  known  properties  of  voltage  and  current 
phasors  In  elementary  electrical  circuit  analysis  that 


2  lZx|2 

av{X<i(t)}  »  -My-l— 


av{Y2(t)>  -  i^L 


av{X(t)Y(t)>  =  - 


Re{ZyZy> 


It  follows  that 


Re{ZxZY) 


p  ■  ’  -  N2/ lz> 


Equation  (94)  Is  completely  equivalent  to  Eq.  (56),  allowing  for  harmless 

scaling;  It  Is  used  In  the  monopulse  radar  literature  to  define  slope 

detection  for  analytical  purposes,  although  It  is  usually  termed  "angle 

detection"  in  recognition  of  the  fact  that  the  slope  measures  the  wave's 

incidence  angle  In  R  units.  However,  the  electronic  circuitry  used  to 

6 

compute  the  angle  Is  often  more  indicative  of  Eq.  (85).  Notice  that  r 
Is  the  slope  of  a  line  drawn  from  the  origin  to  the  point  of  tangency  of 
a  horizontal  line  that  just  touches  the  top  of  the  displayed  ellipse 

A 

(see  Fig.  10).  On  the  other  hand,  R  is  the  slope  of  the  "axis  of 


decorrelation",  l.e. ,  It  Is  the  slope  of  the  ellipse's  major  axis.  To 
the  author's  knowledge,  the  only  justification  for  using  Eq.  (94)  ever 
given  In  the  monopulse  literature  Is  that  Z^/Zy  ought  to  be  real,  and 
that  taking  the  real  part  Is  a  reasonable  way  to  force  a  real  valued 
result  when  the  sinusoids  are  not  perfectly  In  phase. 

It  Is  evident  from  Eq.  (85)  that  r  Is  affected  by  a  scaling  of  the 
coordinates  X  and  Y  In  just  the  same  manner  as  the  slope  of  a  line  would 
be,  l.e.,  It  varies  linearly  with  scaling  of  X  and  Inversely  with  scaling 
of  Y.  The  conversion  of  estimated  slope  to  wave  Incidence  angle  would 
take  any  such  scaling  Into  account  as  a  matter  of  course,  so  the  detected 
wave  angle  would  be  Invariant  to  any  such  scaling  if  r  were  used  to 

A 

estimate  slope.  But  If  R  Is  used  to  estimate  the  slope,  then  the  effect 
of  scaling  Is  somewhat  different.  If  the  X  coordinates  are  scaled  to  a 

A 

very  small  size,  then  p  goes  to  1  and  R  converges  to  r,  and  therefore 
produces  the  same  detected  wave  angle.  On  the  other  hand,  if  the  X 
coordinates  are  scaled  to  larger  values,  the  detected  angle  may  be  sub¬ 
stantially  different,  l.e.,  the  detected  angle  using  R  is  not  invariant 
with  respect  to  display  scaling. 


XIV.  SUMMARY 


The  following  concise  summary  of  processing  methods  Is  presented 
as  an  aid  to  applications.  The  reader  Is  referred  to  the  main  text  for 
details. 


Suppose  Y(t)  Is  driven  by  the  output  of  a  beamformer  having  even- 
symmetric  shading,  and  X ( t )  is  driven  by  the  temporally  integrated  output 
of  a  beamformer  having  odd- symmetric  shading.  If  both  beams  are  formed 
from  the  same  line  array,  then  X(t)  and  Y(t)  are  related  by  a  constant 
multiplicative  factor,  within  certain  restrictions.  For  a  solitary 
sinusoidal  wave  incident  at  angle  e,  they  are  related  by 


where 


X(t)  =  R(e)  Y(t), 

(96) 

K,(e,fo) 

R(e)  =  J2irf  K,(e,f  ) 

(97) 

with  K^(e,fQ)  and  K2(e,fo)  being  the  directivity  functions  of  the  odd  and 


even  beamformers  at  frequency  f 


For  a  solitary  wave  of  arbitrary  spectrum  incident  at  angle  e«l, 
l.e.,  very  near  broadside,  Eq.  (96)  still  holds  with 


R(e)  = 


feV' 


sine 


(98) 


where  ^  is  the  first  moment  of  the  odd  beamformers  aperture  shading 
function  and  M2  Q  Is  the  zeroth  moment  of  the  even  shading  function;  i.e., 
It  Is  the  broadside  sensitivity  of  the  even  beamformer. 
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If  the  even  beamformer  happens  to  be  the  Kerr-Murdock  companion 
to  the  odd  beamformer,  then  Eq.  (96)  holds  for  any  Incident  wave  spectrum 
and  for  all  angles  of  Incidence,  with 

R(e)  *  c"1  sine  .  (99) 

The  locus  of  points  [X(t),Y(t)]  Is  a  straight  line  and  R  Is  Its 
(X  versus  Y)  slope,  1  e.,  the  line  makes  an  angle  3  with  the  vertical 
axis,  where 


6  ■  arctanR(e)  .  (100) 

Over  the  domain  for  which  R(e)  Is  monotonic,  the  angle  of  incidence 
can  be  determined  from  observation  of  R. 

If  several  waves  are  Incident  at  angles  e^Gg,*..,  then  the  locus  Is 
confined  to  the  Interior  of  a  polygon  that  has  two  parallel  sides  for 
each  of  the  slopes  R(e-j),  R(e2),...,  and  the  angle  of  a  component  wave 
can  be  determined  (i.e.,  "resolved")  whenever  the  corresponding  polygon 
sides  are  observed.  In  particular,  for  two  Incident  periodic  waves  the 
locus  will  be  Inscribed  in  a  parallelogram. 

The  predominant  slope  of  a  locus  corrupted  by  noise,  multipath,  or 
other  disturbances  may  be  estimated  as 

R(e) - 2r  ,  (101) 

p  *Vp2  +  4r2 


where 


002) 


(103) 
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2  2 

and  AXy,  Axx,  and  Ayy  are  the  averaged  values  of  XY,  X  ,  and  Y  , 
respectively.  The  average  may  be  computed  by  a  filter,  for  example. 

A 

The  number  r  Is  Itself  an  estimate  of  the  slope,  and  converges  to  R 
If  both  AXy  and  Axx  are  small  compared  with  Ayy. 

A  computational  transform  method  may  also  be  used  to  process  the 
beamformer  outputs,  as  described  In  the  main  text  of  the  paper. 
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APPENDIX 

PERTURBATIONS  OP  APPARENT  DIRECTION 
DUE  TO  PHASE  FRONT  DISTORTION 


be  the  total  field  at  position  r  due  to  spherical,  sinusoidal  waves 
emanating  from  a  collection  of  N  targets.  If  the  targets  are  located  at 
distances  Dj ,Dg,. , . ,DN  from  the  origin,  lying  in  directions  represented 
by  unit  vectors  UpUg,...uN,  then  the  position  of  the  1th  target  Is 
represented  by  the  vector  D^.  Let  W^(r)  denote  the  component  of  the 
complex  envelope  W(r)  due  to  the  1th  target.  Then  it  follows  from  the 
spherical  nature  of  each  wave  tnat 

Wl(r).  V0)(^-)e*Pp^  .  (A.2) 

where  D.  denotes  the  distance  between  the  point  r  and  the  ith  target, 
ir 

D.  &  DjU.-r  .  The  field  in  the  vicinity  of  the  origin  may  be  found  by 

i  r  ii 

expanding  Eq.  (A.2)  to  first  order  in  the  differential  vector  dr, 

(dr)  -  W.(0)(l+jK1<ui,dr>)  ,  (A. 3) 

where 

2irf  , 

K.  A  -_o  .  j-  (A. 4) 

i  c  D1 

and  <a,b>  denotes  the  inner  product.  Notice  that  quickly  converges 
to  the  real  valued  wave  number  constant,  2iTf0/c,  as  the  target  distance 
Increases.  The  total  field  due  to  the  N  targets  is 

W(dr)  =  £  W.(dr)  =  (l+j<G,dr>)W(0)  ,  (A. 5) 

i  1 
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‘WtllHflWiP: 


where 


z  Wi(0)K1u1 

G  -  1  z  W, (0)  *  (A,6) 

i  1 

and  we  have  made  use  of  the  fact  that  W(0)BW.j (0) .  Contemplation  of 
Eq.  (A. 5)  leads  to  the  conclusion  that  Re{6)  is  the  gradient  vector  of 
the  phase  of  the  field  in  the  vicinity  of  the  origin,  while  -|W(0)|Im{G} 
Is  the  gradient  vector  of  the  amplitude. 

Let  us  consider  some  special  cases.  Suppose  there  is  just  one 
target  at  distance  D-j  in  direction  u^.  Then  G=K1.u^,  so  the  phase 
gradient  is  (2irfQ/c)ui.  while  the  amplitude  gradient  is  | W(0)  | u-j /D^ . 

Both  gradient  vectors  point  in  the  direction  of  the  target,  i.e.,  towards 
the  source  of  the  single  spherical  wave.  The  amplitude  gradient  is  very 
weak  for  a  distant  target,  but  the  phase  gradient  remains  strong  even  at 
extreme  distances. 

The  outcome  is  similar  if  there  are  N  targets  lying  in  the  same 

direction,  u^u^ . uN,  since  G  is  then  a  scalar  multiple  of  u] .  Both 

gradient  vectors  point  in  that  direction  and  the  magnitude  of  the  phase 
gradient  vector  still  approximates  the  wave  number  Zvfjc,  provided  all 
of  the  targets  are  rather  distant. 

Now  suppose  there  are  N  targets  at  different  directions  but  with 
positions  and  phase  characteristics  that  connive  to  make  the  spherical 
waves  arrive  at  the  origin  in  phase.  Their  shared  phase  component  can 
then  be  canceled  in  the  numerator  and  denominator  of  Eq.  (A. 6).  In  this 
case,  the  phase  gradient  is 


27rfQ  ?  Iwi(°)lui 
Re{G}  =  — -  j;  |  wi  (0)  | 


(A.  7) 
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This  formula  Is  convex  with  respect  to  the  u^'s  because  the  weighting 
coefficients  are  all  positive;  i.e.,  the  phase  gradient  lies  within  a 
convex  subset  bounded  by  the  vectors  u^,  and  points  to  a  particular  spot 
somewhere  within  the  collection  of  targets,  termed  the  "radar  centroid" 
in  monopulse  literature.  The  amplitude  gradient  vector  has  a  similar 
form: 


z  |Wi|ui/D1 

-|w(o)|  im (G >  =  | w(o) |  —  *  (A,8) 

If  all  the  targets  are  at  roughly  the  same  distance  D  ,  then  the  ampli- 

L> 

tude  gradient  is  approximately  (c/27rf_) |W(0) | D  times  the  phase  gradient, 

0  V 

so  they  point  in  virtually  the  same  direction. 

When  the  spherical  waves  do  not  arrive  in  phase,  the  outcome  is  much 
more  complicated.  For  active  sonar  applications,  we  are  most  interested 
in  the  case  of  a  fairly  distant  collection  of  N  targets  that  are 
clustered  around  an  arbitrarily  defined  geometric  center  whose  distance 
and  direction  from  the  origin  are  represented  by  Dc  and  the  unit  vector 
u  ,  respectively.  We  assume  the  cluster  has  a  limited  downrange  extent 

V 

in  the  sense  that  D^/Dc  is  extremely  close  to  unity  for  all  i,  although 
the  differences  are  permitted  to  be  much  larger  than  the  acoustic 

wavelength.  In  effect,  this  means  (2rrf  /c-j/D  )  can  be  substituted  for 

U  V 

the  K.|'s  in  Eq.  (A. 6),  with  the  result  that 


s  (X-.+jYJu,  *  K  £ 
i  1  1  1  c  i 


Ziui 


(A. 9) 


to  a  very  close  approximation,  where 


and 


2lff  -i 

K_  4  — -  4- 


(A.ll) 


z  (0)  so  z  X.  s  1  and  z 
J^crossrange^extent  then  ih 


Y1  -  0. 


If  the  target 


e  direction  vectors  will 


Notice  that  W(0) 
cluster  has  limited 'crossrange 

all  lie  very  close  to  ur,  and  it  follows  that  Re{U)  =  z  X4u„  =  1 ,  and 

•W  ^  ^  I  C 

Im{U}  «  z  Y.u  =  0,  i.e.,  U  =  u  .  We  can  investigate  the  variation  of  LI 

■j  '  V  t 

from  thii  approximate  value  by  setting  ^=uc+(Uj-uc)  in  Eg.  (A. 10), 


u  •  uc  +  U 


(A. 12) 


where  the  complex  displacement  vector  U  is  defined  as 


U  4  r  Zi(ui-uc)  .  (A. 13) 

Although  uc  and  the  u^'s  are  real,  and  Kc  is  almost  purely  real,  there 
is  no  simple  rule  for  predicting  the  phase  of  the  displacement  vector  U 
because  the  complex  ratios  can  have  arbitrary  phases.  Even  the 
direction  of  U  is  uncertain,  since  it  is  a  linear  combination  of  the 
vectors  u^-uc,  which  are  likely  to  point  in  a  variety  of  directions  if 
the  targets  are  clustered  randomly  about  the  geometric  center. 

The  amplitude  gradient  vector  points  in  the  direction  of  -Im(G), 
which  can  be  expressed  using  Eqs.  (A. 10)  and  (A. 12)  as 

u. 

-Im(G)  =  -Im{KcU)  =  ^  -  Im{KcU>  .  (A. 14) 
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Its  direction  Is  unpredictable  due  to  the  aforementioned  uncertainty  in 
U,  which  is  likely  to  be  the  dominant  factor  since  the  term  uc/Dc  is  so 
small  for  a  distant  target.  There  is  little  that  can  be  said  about  the 
amplitude  gradient  in  this  case,  other  than  that  It  is  likely  to  be  very 
weak  and  may  not  point  in  the  same  direction  as  the  phase  gradient. 
Fortunately,  the  amplitude  gradient  is  an  unimportant  factor  here.  Based 
upon  radar  experience,  any  odd-even  beamformer  pair,  any  monopulse  radar, 
or  virtually  any  tracking  radar  will  characteristically  perceive  the 
target  cluster  to  be  in  the  direction  of  the  phase  gradient  vector. 

To  get  the  phase  gradient  of  the  acoustic  field  we  must  take  the 
real  part  of  G, 


Re{G}=  Re{K  G}  = 

v 


(uc+Re{U>)  + 


(A. 15) 


The  last  term  can  be  neglected  since  the  target  cluster  is  distant,  so 
the  phase  gradient  points  in  the  direction  of  Re{U}=uc+Re{U}.  We  may 
think  of  it  as  pointing  toward  the  apparent  center  of  the  target  cluster, 
also  termed  the  phase  center  in  the  radar  literature.  The  apparent  cen¬ 
ter  is  displaced  from  the  geometric  center  in  accordance  with  Eq.  (A. 15). 
The  quantities  in  the  definition  of  U  given  by  Eq.  (A. 13)  express  the 
fractional  contribution  of  each  target  to  the  total  field  at  the  origin 
W(0),  and  can  be  expected  to  be  of  less  than  unit  magnitude  under  normal 
circumstances.  Due  to  the  limited  crossrange  extent  of  the  cluster,  the 
vectors  (u--ur)  are  extremely  small.  Thus  the  vector  U  must  be  of  very 
small  magnitude,  and  the  apparent  center  is  only  slightly  displaced  from 
the  geometric  center.  However,  if  the  configuration  and  orientation  of 
the  cluster  are  such  as  to  produce  destructive  interference  of  the  com¬ 
ponent  spherical  waves  at  the  origin,  i.e.,  if  the  origin  lies  in  the 
direction  of  a  "null"  in  the  radiation  pattern  of  the  cluster,  then  W(0) 
may  be  small,  even  zero.  Under  that  circumstance,  the  fractions 
Zi-Wi(0)/H(0)  may  become  very  large,  making  the  displacement  vector  U 
quite  significant.  Indeed,  the  apparent  center  may  jump  entirely  outside 
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the  target  cluster,  and  its  location  may  fluctuate  erratically  when  the 
cluster's  orientation  or  configuration  Is  changed.  In  early  radar 
experience,  this  effect  was  often  encountered  when  an  aircraft  target, 
which  can  be  regarded  as  a  cluster  of  point  scatterers,  performed  rapid 
maneuvers  or  reacted  to  turbulent  air.  The  effect  may  have  first  been 
attributed  to  brief  flashes  of  specular  radar  reflection  at  different 
positions  on  the  aircraft  caused  by  rotation  thus  explaining  the  use  of 
the  term  glint.  Theoretical  analyses  refute  this  explanation,  however, 
revealing  that  the  effect  actually  occurs  »iran  the  received  wave  Is 
weakest. 

In  summary,  for  the  case  of  a  distant  cluster  of  N  targets,  the 
phase  gradient  of  the  total  field  at  the  origin  points  in  the  direction 
of  the  vector  \ 


Re{U}  =£  X,u.  ,  (A. 16) 

i  1  1 

where  u^  Is  the  unit  vector  pointing  toward  the  ith  component  target  and 
the  X.j  Is  the  real  part  of  the  complex  ratio  of  Its  contributed  field  to 
the  total  field  at  the  origin,  with  the  X^'s  obeying  i  X^  =  1.  The 
vector  Re{U)  defines  the  apparent  direction  of  the  target  cluster,  1 .e . , 
the  direction  normal  to  the  incident  phase  front.  It  Is  clear  from 
Eq.  (A. 16)  that  the  vector  Re{U}  will  be  attracted  toward  any  exception¬ 
ally  strong  component  of  the  cluster,  and  will  at  least  point  toward  the 
interior  of  the  cluster  unless  there  Is  enough  destructive  Interference 
of  the  waves  that  the  glint  phenomenon  intrudes. 

Consider  which  apparent  directions  are  theoretically  possible  for  a 
target  cluster  of  a  given  geometry  if  we  permit  the  phase  and  amplitude 
contributed  by  each  component  to  be  varied  independently.  We  observe 
that,  given  any  numbers  a^  that  have  a  positive  sum  S,  we  could  theo¬ 
retically  adjust  the  phase  and  amplitude  of  each  target's  wave  to  give 
W.j (0)=a1* .  Then  W(0)  *  £  W^O)  =  S,  and  X^a^Oj/S.  The  apparent 
direction  Is  that  of  the  vector  RefU}®  S"^  £  a^.  We  thus  conclude  that 
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the  apparent  direction  can  correspond  to  any  linear  combination  of  the 
u^'s,  provided  the  weighting  coefficients  have  a  positive  sum.  This  con 
elusion  Is  almost,  but  not  quite,  equivalent  to  saying  that  the  apparent 
direction  can  lie  anywhere  In  the  subspace  spanned  by  the  u^'s.  For  a 
two-target  cluster,  the  apparent  direction  can  be  anywhere  in  a  half¬ 
plane  containing  the  cluster  whereas  for  a  three-target  cluster  It  can 
be  anywhere  in  a  half-volume  which  contains  it.  As  an  example,  consider 
a  two-target  cluster.  Although  It  may  seem  that  there  are  four  Indepen¬ 
dent  parameters  to  deal  with,  namely  the  real  and  Imaginary  parts  of 
W-j (0)  and  W2(0),  it  is  really  only  necessary  to  consider  =Re{W1  (0) > . 
The  phase  gradient  is  27rf0/cxRe{U},  and  Re{U}=X^u^X2U2=U2+X-j (u-|-u2) , 
since  X-j+X2=l.  So  the  apparent  direction  Is  that  of  a  vector  whose  end 
point  lies  on  the  straight  line  passing  through  the  end  points  of  the 
unit  vectors  u^  and  u2  that  point  toward  the  two  component  targets.  If 
X-|  and  X2  are  both  positive,  then  that  vector  lies  between  u^  and  u2» 
but  if  the  waves  destructively  Interfere,  i.e.,  if  X-|  and  X2  have 
opposite  signs,  then  the  vector's  end  point  may  lie  elsewhere  along  the 
line.  Notice  that  the  phase  gradient  is  of  largest  magnitude  when  its 
deviation  is  greatest,  and  when  the  destructive  Interference  of  the 
waves  is  most  severe.  Although  this  example  Illustrates  that  extreme 
deviations  in  apparent  directions  of  a  target  cluster  are  possible,  they 
are  not  necessarily  probable.  The  use  of  the  parallelogram  sketching 
method  as  described  in  the  main  text  of  this  report  may  be  of  value  in 
predicting  such  effects. 

The  perturbations  in  the  direction  of  the  phase  gradient  are 
symptoms  of  "kinks"  in  the  phase  fronts  that  radiate  from  the  cluster. 
Since  they  occur  at  points  where  destructive  interference  takes  place, 
they  must  be  aligned  with  the  "null"  directions  in  the  directivity  pat¬ 
tern  of  the  cluster.  The  target  cluster  may  be  regarded  as  a  trans¬ 
mitting  array  having  an  unusual  shape.  Figure  A.l  shows  phase  fronts 
emanating  from  a  two- target  cluster  with  three-wavelength  spacing,  in 
which  the  emanated  waves  are  of  almost  equal  strength  but  are  driven 
180°  out  of  phase. ^  Significant  perturbations  in  the  wave  fronts  are 
seen  at  only  eight  specific  aspect  angles.  (Destructive  interference 


FIGURE  A.1 

DISTORTED  PHASE  FRONTS  FROM  TWO  WAVE  CENTERS  SEPARATED 
BY  THREE  WAVELENGTHS  AND  DRIVEN  180°  OUT  OF  PHASE 
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also  takes  place  at  the  end-fire  angles,  but  without  any  attendant 
distortion  In  the  phase  front.)  Note,  however,  that  In  the  case  of 
active  sonar  or  radar  Illumination  from  a  fixed  direction,  a  simple 
rotation  of  the  cluster  does  not  produce  an  Identical  rotation  of  the 
kink  pattern.  The  actual  result  Is  made  more  complicated  by  the  changing 
phase  relationship  of  the  field  Incident  upon  the  target  components. 
Extensive  measurements  of  the  glint  properties  of  specific  cluster  geom¬ 
etries  (e.g.,  aircraft)  have  been  performed  in  radar  laboratories. 

The  treatment  given  in  this  appendix  has  used  vectors  rather  than 
coordinate  angles  to  represent  target  directions,  but  it  may  readily  be 
verified  that  the  bearing  angle,  relative  to  any  "horizon,"  of  the  vector 
iX.jUj,  where  the  u^'s  are  almost  collinear  unit  vectors  and  the  numbers 
L  sum  to  unity,  is  given  to  first  order  by  z  X.e.,  where  is  the  bearing 

1  Jill 

angle  of  u^.  For  the  case  of  interest,  X1-=Re{W^ (0)/W(0) }  so  the  apparent 
bearing  angle  of  the  target  cluster  ec  is  given  to  first  order  by 


e 


c 


Re 


W^O) 


wW 


9.j  =  Re 


£0^(0) 

Iw.(o) 

i  1 


(A. 17) 


This  expression  is  almost  the  same  as  that  given  in  Eq.  (79)  for  the 
sensing  ratio  used  in  monopulse  radar,  repeated  here  for  convenience. 


R  =  ^ -  .  (A. 18) 

Lq?.i 

The  two  formulas  are,  in  fact,  equivalent  for  a  small  cluster  near 
broadside.  Recall  that  the  R's  measure  the  bearing  angles  in  R  units, 
which,  for  targets  near  broadside,  give  a  linear  measure  of  angle,  and 
the  constant  of  proportionality  cancels  in  numerator  and  denominator. 
Furthermore,  if  the  cluster  is  small  with  respect  to  the  undulations  of 
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the  beam  pattern  of  the  even  beamformer  of  an  odd-even  pair,  then 
Q2  K2(ec)W1(0),  where  K2(ec)  Is  the  even  directivity  pattern  function 
evaluated  In  the  direction  of  the  cluster,  and  this  factor  also  cancels 
In  numerator  and  denominator.  The  conclusion  Is  that  an  odd-even  beam- 
former  pair  perceives  the  cluster  to  lie  In  the  direction  of  the  phase 
gradient,  perpendicular  to  the  Incoming  phase  front.  The  same  conclu¬ 
sion  applies  to  monopulse  radar,  as  has  been  noted  In  the  literature.'®*'5 


86 


REFERENCES 


1.  W.  6.  Bruner  et  al . ,  "Radar  and  Navigational  Systems,"  Electronic 
Designers  Handbook,  edited  by  L.  J.  Giacoletto  (McGraw-Hill  Book  Co., 
Inc.,  New  York,  1977) . 

2.  P.  Z.  Peebles  and  R.  S.  Berkowltz,  "Multiple-Target  Monopulse  Radar 
Processing  Techniques,"  IEEE  Trans.  Aerospace  &  Electronics  Systems, 
845-854,  (1968). 

3.  J.  Capon,  "High-Resolution  Frequency-Wavenumber  Spectrum  Analysis," 

Proc.  IEEE,  1408-1418  (1969). 

4.  J.  S.  Kerr  and  W.  L.  Murdock,  "Relations  Between  the  Far  Field  and 
the  Illumination  of  Radar  Apertures,"  GE  Report  No.  51-E-234, 

General  Electric  Company,  1  November  1951,  Appendix  B. 

5.  G.  M.  Kilpatrick,  "Final  Engineering  Report  on  Angular  Accuracy 
Improvement,"  General  Electric  Company  Report,  1  August  1952, 

Appendix  I.  (See  Ref.  9.) 

6.  D.  R.  Rhodes,  Introduction  to  Monopulse  (McGraw-Hill  Book  Co.,  Inc., 

New  York,  1959T 

7.  Special  Issue  on  Time  Delay  Estimation,  IEEE  Trans.  Acoustics,  Speech, 

&  Signal  Processing,  July  1981. 

8.  National  Defense  Research  Committee,  "The  Bearing  Deviation  Indicator" 
in  "Underwater  Sound  Equipment  II,  Echo-Ranging  Systems,"  Summary 
Technical  Report  of  Division  6,  NDRC,  Office  of  Scientific  Research  and 
Development,  Washington,  D.C.,  1946,  Vol.  15,  pp.  83-116. 

9.  D.  K.  Barton  (ed.),  Radars,  Vol.  1,  Monopulse  Radar  (Artech  House,  Inc., 
1974). 

10.  R.  A.  Ross  and  M.  E.  Bechtel,  "Scattering  Center  Theory  and  Radar 
Glint  Analysis,"  IEEE  Trans.  Aerospace  &  Electronic  Systems,  756-762 
(1968). 

11.  I.  Kanter,  "Varieties  of  Average  Monopulse  Responses  to  Multiple 
Targets,"  IEEE  Trans.  Aerospace  &  Electronic  Systems,  25-28  (1981). 

12.  D.  K.  Barton,  "Accuracy  of  Monopulse  Radar,"  IRE  Proceedings  of  the 
3rd  National  Convention  on  Military  Electronics,  1959.  (See  Ref.  9.) 


87 


13.  J.  W.  Horton,  Fundamentals  of  Sonar  (United  States  Naval  Institute, 
Annapolis,  Md.,  1959). 

14.  S.  M.  Sherman,  "Complex  Indicated  Angles  Applied  to  Unresol vable 
Radar  Targets  and  Multipath,"  IEEE  Trans.  Aerospace  &  Electronic 
Systems,  160-170,  1971. 

15.  D.  0.  Howard,  "Radar  Target  Angular  Scintillation  In  Tracking  and 
Guidance  Systems  Based  on  Echo  Signal  Phase  Front  Distortion," 
Proceedings  of  National  Electronics  Conference,  Vol.  15,  1959. 

16.  M.  G.  Kendall  and  A.  Stuart,  Chapter  29  of  The  Advanced  Theory  of 
Statistics  (Charles  Griffin  &  Co.,  Ltd.,  London). 


88 


17  Hay  1982 


DISTRIBUTION  LIST  FOR 
ARL-TM-82-10 

UNDER  CONTRACT  N00024-79-C-6358,  Task  8,  Subtask  14 

UNCLASSIFIED 


Copy  No. 


1 

2 

3 

4 

5 

6 

7 

8 
9 

10 


11 

12 

13 

14 

15 


16 


17 

18 


19 


Commander 

Naval  Sea  Systems  Command 
Department  of  the  Navy 
Washington,  DC  20362 
Attn:  '  0.  Porter  (Code  63R1) 

A.  Carpenter  (PMS  407) 

C.  Smith  (Code  63R) 

E.  Liszka  (PMS  406) 

O.  Neely  (PMS  409C) 

F.  Romano  (Code  63R3) 

B.  Chapman  (Code  63Z34) 

L.  Epperly  (Code  63D) 

R.  Snuggs  (PMS  409) 

Library  (Code  9961) 

Commanding  Officer 

Office  of  Naval  Research 

Arlington,  VA  22217 

Attn:  A.  Ellenthorp  (Code  220) 

R.  Obrochta  (Code  230) 

CAPT  A.  Gilmore  (Code  220) 

P.  Rodgers  (Code  425AC) 

L.  Hargrove  (Code  412) 

Commanding  Officer 

Naval  Research  Laboratory 

Underwater  Sound  Reference  Detachment 

P.  0.  Box  8337 

Orlando,  FL  32856 

Attn:  J.  Blue 

Commanding  Officer 

Naval  Ocean  Research  and  Development  Activity 
NSTL  Station,  MS  39529 
Attn:  E.  Chaika 

B.  Blumenthal 

Defense  Advanced  Research  Projects  Agency 
1400  Wilson  Boulevard 
Arlington,  VA  22209 
Attn:  T.  Kooij,  TTO 


89 


Distribution  List  for  ARL-TM-82-10  under  Cont.  N00024-79-C-6358,  Task  8, 
Subtask  14 

Copy  No. 

Naval  Electronic  Systems  Command 
Department  of  the  Navy 
Washington,  DC  20360 

20  Attn:  R.  Mltnlk,  Code  612 

21  Commander 

Naval  Material  Command 
Department  of  the  Navy 
Washington,  DC  20360 

Commander 

Naval  Ocean  Systems  Center 
Department  of  the  Navy 
San  Diego,  CA  92152 

22  Attn:  J.  Reeves 

23  R.  Anderson 

24  R.  McLennan 

25  H.  Sexton 

Commanding  Officer 

Naval  Oceanographic  Office 

NSTL  Station,  Bay  St.  Louis,  MS  39522 

26  Attn:  W.  Jobst 

27  G.  Lewis 

Commanding  Officer 

Naval  Coastal  Systems  Center 

Panama  City,  FL  32407 

28  Attn:  L.  Flax 

29  S.  Richardson 

Commander 

Naval  Surface  Weapons  Center 
White  Oak  Laboratory 
Silver  Spring,  MD  20910 

30  Attn:  Library 

Superintendent 

Naval  Postgraduate  School 

Monterey,  CA  93940 

31  Attn:  Library 

32  Commanding  Officer  and  Director 
Defense  Technical  Information  Center 
Cameron  Station,  Building  5 

5010  Duke  Street 
Alexandria,  VA  22314 


90 


Distribution  List  for  ARL-TM-82-10  under  Cont.  N00024-79-C-6358,  Task  8, 
Subtask  14 


33 

Applied  Physics  Laboratory 

University  of  Washington 

1013  N,E.  40th  Street 

Seattle,  WA  ru!05 

Attn;  L1br*"v 

34 

Applied  Research  Laboratory 

The  Pennsylvania  State  University 

P.  0.  Box  30 

State  College,  PA  16801 

Attn;  Library 

35 

3b 

37 

Admiralty  Underwater  Weapons  Establishment 
Portland,  Dorset  DT5  2 JJ 

UNITED  KINGDOM 

Attn;  R,  Bugler 

0.  Weston 

J.  Willis 

38 

39 

University  of  Rhode  Island 

Kingston,  RI  02881 

Attn:  Prof.  D.  W.  Tufts 

Prof.  L.  Scharf 

40 

Yale  University 

New  Haven,  CT  06520 

Attn:  Prof.  P.  M.  Schulthelss 

41 

Princeton  University 

Princeton,  NO  08540 

Attn:  Prof.  S.  Schwartz 

42 

University  of  Mississippi 

University,  MS  38677 

Attn:  Prof.  L.  W.  Pearson 

43 

Electronic  and  Electric  Engineering  Division 
1'ilverslty  of  Technology 

Loughborough,  Leicestershire 

UNI' .0  KINGDOM 

Attn;  R.  Griffiths 

44 

45 

Department  of  Physics 

The  University  of  Texas  at  Austin 

Austin,  TX  78712 

Attn:  W.  Nolle 

T.  Griffy 

91 


Distribution  List  for  ARL-TM-82-10  under  Cont.  N00024-79-C-6358,  Task  8, 
Subtask  14 

ggREA, 


46 

47 

48 

49 

Department  of  Electrical  Engineering 
The  University  of  Texas  at  Austin 
Austin,  TX  l2 

Attn:  E.  P<  s 

F .  Bo  iwk 

T.  Itoh 

E.  Hixson 

50 

Office  of  Naval  '«rc^ 

Resident  Represen  Iv 

Room  No.  582 

Federal  Building 

Austin,  TX  78701 

51 

Electroacoustics  Group,  ARL:IIT 

52 

Signal  Physics  Group,  ARL:UT 

53 

Environmental  Sciences  Group,  ARL:UT 

54 

Sonar  Engineering  Group,  ARL:UT 

55 

Physical  Sciences  Group,  ARL :UT 

56 

Sonar  Development  Division,  ARL : UT 

57 

Lewie  M.  Barber,  ARL:UT 

58 

Garland  R,  Barnard,  ARL :UT 

59 

George  P.  Coble,  ARL : UT 

60 

Marshall  E.  Frazer,  ARLilIT 

61 

Terry  L.  Henderson,  ARL : UT 

62 

Robert  E.  Hollingsworth,  ARL :UT 

63 

Arthur  C.  Holly,  ARL:UT 

64 

Claude  W.  Horton,  ARL:UT 

65 

John  M.  Huckabay,  ARL : UT 

66 

John  W.  Maxwell,  ARL :UT 

67 

Chester  M,  McKinney,  ARL : UT 

68 

T.  G,  Muir,  ARL ;UT 

92 


Distribution  List  for  ARL-TM-82-10  under  Cont.  N00024-79-C-6358,  Task  8 
Subtask  14 

Copy  No. 

69  S.  Patrick  Pitt,  ARL-.UT 

70  Dennis  R.  Powell,  ARL:UT 

71  Robert  H.  Stokes,  ARL:UT 

72  Reuben  H.  Wallace,  ARL : UT 

73  Joseph  F,  Willman,  ARL :UT 

74  Gary  R.  Wilson,  ARL :UT 

75  Library,  ARL:UT 
SEG  Reserve,  ARL :UT 


76  -  96 


